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A Area of a structural member 
Tributary area of a structur- 
(sq. ft) 

A_ Area of steel beam cross sec 
(sq. in.) 

A. Net area of the rear face of 
^ (sq. ft) 

A Area of concrete in cross se 
c 

A- Area of flange of structural 
Ifet area of the exterior fro 
vith openings (sq. ft) 

A Cross section area of concre 

g 

A - Gross area of tne exterior f 

* (sq. ft) 

A Area of each portion of the 
n ture vith openings (sq. ft) 

A - Area of openings in the f ron 
01 / .\ 
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a 1 Distance from edge "a" to < 
a' T Distance from edge "a" to 

B Peak value of externally a] 
Elastic stiffness coefficit 

B Equivalent peak value of tl 

B Peak value of the effective 
^ multi-story building 

B Peak value of the external 

6 

b Width of flange of steel be 
Width of compression flange 
Long side dimension of two- 

b Flange width of channel sec 
b Width of longitudinal web 

b 1 Width of stem of reinforced 
Width of flange of steel b 
Distance from edge "b" to c 

b" Distance from edge "b" to t 
C Compression force developed 

n rrrnvnvp' .Q n 1 on rnlnmn i 



Symbols 



c Distance from neutral axis to 
Cohesive strength of the soil 
laboratory testing procedures 
The clear distance betveen coi 

c Velocity of sound in undisturl 
c _ Velocity of sound in the regi< 

D Diameter of a cylindrical segj 
Diameter of the smallest long: 
Diameter of spiral core of a i 

D.L.F Dynamic Load Factor = x /x 
J m' s 

D.I.F Dynamic Increase Factor for s* 

d Depth of burst of bomb below j 
Height of burst of bomb above 
Depth from the compression fa< 
of longitudinal tensile reinfi 
Least lateral dimension of re< 
Total depth of steel beam (in 
Over-all depth of truss (ft) 
Diameter or width of column cs 
Diameter of circular reinforci 

d' Distance betveen centroids of 

Hmi'hlv Tf^i TTPnvr^nl rn < nrT*a'f.#a moi 



EM 

15 *fer 57 

F Summation of all external horizoi 

including foundation reactions 

F Pressure factor for underground 1 

^ Normal component of total passiv< 

F(t) Load applied to a structural ele'i 

f Fiber stress (psi) 

f Column critical "buckling stress 

1 

f , Average stress at design level di 
ad 

f Dynamic yield strength for struc 
construction only (psi) 

f Elastic buckling stress due to b< 

f,. Maximum stress due to bending at 
subjected to axial load and bend 

f Dynamic yield strength of steel 
ay . , 

f (t) Effective load on the g floor 

^ function of time 

f. Initial static stress (psi) 

f. Lover static stress (psi) 

P TVimftTn-1 f* xHtfaTH c+Tpmcrh)") for* T*^"! nf* 



Symbols 



H Equivalent inrpulse acting on tl 

H , Total impulse of the external : 
gt 

H, Total impulse of the external : 

v 

h Height (ft) 

Story height (ft) 

Column height (ft) 

Height of a structure above gr< 

Maximum thickness of channel f! 

h Clear height of column (ft) 
c . 

h Height of the g story of a nn 

o 

h Clearing dimension for each poi 
of rectangular structure vith < 

h f Clearing height of a closed re< 
Vertical distance from roof to 
structure (ft) 

h} Weighted average buildup heigh" 
rectangular structure with opei 

h' Weighted average clearing heig] 
a rectangular structure vith O] 

* Weighted average 
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KE 

(KE) a 

(KE) 



KLM 



s 

KT 

K 
K 1 

K" 
k 



ep 

*!' 



g 



Kinetic energy- 

Kinetic energy of the actual system 

Kinetic energy of the equivalent system 

Either load, mass, or load-mass factor for slab fixed on four sides 

(i = 1, 2, 3,...; j = 1, 2, 3 ---- ) Stiffness influence coefficient 
(kips/ft) 

Load factor 
Load-mass factor 
Mass factor 

Normal component of passive pressure coefficient accounting for 
the cohesive effect of any soil having an internal friction equal 
to 

Normal component of passive pressure coefficient for any soil with 
internal friction angle equal to and with zero wall friction 
developed 

Normal component of passive pressure coefficient for any soil with 
internal friction angle equal to and with maximum wall friction 
developed 

Resistance factor 

Either load, mass, or load-mass factor for slab simply supported 
on four sides 

Kiloton, 1000 tons 

Either load, mass, or load-mass factor for special edge conditions 

Ratio of the maximum average overpressure to the reflected over- 
pressure existing on an inclined roof 

Beam equivalent length coefficient used in the design of steel 
beam columns 

Column length factor used in the design of steel columns and beam 
columns 

Spring constant, force required to cause unit deflection of spring 

(Hips/ft) 

Soil pressure factor (psi) 

3oil constant for underground explosion (psi) 
Effective spring constant (kips/ft) 
Equivalent spring constant (kips/ft) 
Gpring constant in the elasto-plastic range 



*i t k , 



...k g ... 1'prins constants for the first, second, and the g th 
stories of a multi-story building (kips/ft) 



Simplified notation of k 



g,g-l 
X 



, 
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k Spring constant of the coupling spring between the g and the i th 
6 floors of a multi-story building (kips/ft) 

k', k' = - k . 
gi gi gi 

k cl' k K2'" k Ki'"' k sm s P rin S constants of the springs connecting m with 

the masses m.., m_, . . .m,, . . .m g 

L Span length of beam or truss (ft) 
Unsupported length of beam (ft) 

Length of a rectangular structure in the direction of propagation 
of the blast wave (ft) 

Length of a cylindrical segment along the cylinder axis (ft) 
Length of shear wall center to center of column steel 
Spacing of columns in each direction 

L. Distance from the outside of the front face to the inside rear 
face of a rectangular structure with openings (ft) 

L 1 Distance from the front face of a rectangular structure to a point 
under consideration on the roof or sides in the direction of pro- 
pagation of the blast wave (ft) 

Length between sections of zero and maximum moment being considered 
(composite beams) 

M Bending moment applied to a section 

Moment of forces on piles about their centroidal axis 
Resisting moment of soil on the footing per unit width 

M Bending moment at center line of a beam or slab 

M Negative resisting moment in column per foot, flat slabs (kip-ft) 

M Positive resisting moment in column strip per foot, flat slabs 
^ (kip-ft) 

M Maximum design moment in a member under axial load, P 

M, Dynamic bending moment 

M , Elastic dynamic buckling moment 

M Maximum bending moment 
m 

M- Component in a plane perpendicular to edge "a" of the total re- 
sisting moment along the fracture lines bounding area "A," two- 
way slabs (kip-ft) 

M-. Component in a plane perpendicular to edge "b," of the total re- 
sisting moment along the fracture lines bounding area B, two-way 
slabs (kip-ft) 

M Negative resisting moment in middle strip per foot, flat slabs 
11111 (kip-ft) 

M Positive resisting moment in middle strip per foot, flat slabs 
** (kip-ft) 

XI 
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Negative "bending moment at support 

Moment of all external forces on structure about axis of rotation 



Pcn 



"pfa 



"pfb 



j?np 

M 
pos 

"PS 

"Psa 
"Psb 



"pmb 



Summation of external moments about the point of rotation ex- 
cluding footing -projection moments 

Plastic resisting moment under "bending only 

Negative plastic resisting moment in column strip per foot, flat 
slabs (kip-ft) 

Positive plastic resisting moment in column strip per foot, flat 
slabs (kip-ft) 

Component of the total plastic bending moment capacity along the 
fracture line boundary oT area "A" which is in a plane perpendic- 
ular to edge "a," or the total positive plastic bending moment 
capacity for a section parallel to edge "a," two-way slabs (kip-ft) 

Component of the total plastic bending moment capacity along the 
fracture line boundary of area "B" which is in a plane perpendic- 
ular to edge "b," or the total positive plastic bending moment 
capacity for a section parallel to edge "b," two-way slabs (kip/ft) 

Plastic resisting moment at centerline of beam or slab 

Negative plastic resisting moment in middle strip per foot, flat 
slabs (kip-ft) 

Positive plastic resisting moment in middle strip per foot, flat 
slabs (kip-ft) 

Maximum positive bending moment (kip-ft) 
Plastic resisting moment at support (kip-ft) 

Total negative plastic bending' moment capacity along edge "a, " 
two-way slabs (kip-ft) 

Total negative plastic bending moment capacity along edge "b," 
two-way slabs (kip-ft) 

Yield resisting moment (kip-ft) 

Plastic resisting moment per unit of width of slab (kip-ft/ft) 

Plastic positive bending moment capacity per unit width for short 
span, two-way slabs (kip-ft/ft) 

Plastic negative bending moment capacity per unit width at center 
of edge a for long span, two-way slabs (kip-ft/ft) 

Plastic negative bending moment capacity per unit width at center 
of edge b for short span, two-way slabs (kip-ft/ft) 
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M, Moment at the intersection of the two-linear portions of the P vs ! 
curve of a steel beam section (kip-ft) 

Ml Theoretical plastic resisting moment (kip-ft) 

{ 2 n \ 
m Mass per unit length ^kip-sec /f t c ) 

Point mass ^kip-sec /ft) 

In concrete design: f /0.85 f^ c 

Total moving mass of structure and earth included between footings 
Number of fundamental dimensional quantities 

m Mass of the equivalent system y kip-sec /ft ) 

m Mass of structure considered to rotate as well as translate 

m. Total mass of the element or structural system under consideration 

(kip-sec 2 /") 
m. Total mass of the equivalent system (kip-sec /ft) 

TV\ 

m.. , nu. .m Mass of the first, second, and g floors of a multi-story 
1 d g building 

N Number of non-dimensional parameters 

Total number of shear connectors required between the points of 
zero and maximum moment of a composite beam 
Weighted number of piles in group 

N_ Ratio of the length of fixed-edge perimeter to the total perimeter 

N Ratio of the length of simply -supported perimeter to the total 
perimeter 

n Number of columns in a story 

Ratio of modulus of elasticity of steel to modulus of elasticity 

of concrete 

Number of dimensional variables 

Number of stories in a multi-story frame 

P Load or force (kips) 

Total dynamic load on slab (kips) 
Panel influence factor for shear walls 

P. Total load on "A" portion of two-way slab (kips) 

r\ 

P Average axial load acting on each of several columns of a frame 
&V (kips) 

P Local overpressure on the back face of a buried rectangular 
structure (psi) 

P_ Total load on "B" portion of two-way slab (kips) 
B 

P Compression mode overpressure (psi) 

P Uniform compressive pressure applied radially on arch 

Cl 

P Uniform radial pressure that produces buckling of arch 

XIII 



EM 
15 Mar 

P cyl 



Symbols 



57 



dome 



*ed 

P e (t) 

? f 
P front 



Local overpressure normal to the exterior surface of a cylindrical 

segment (psi) 

Maximum axial load on column with given M^ (kips) 

Reduced maximum axial load for long columns 

Deflection mode overpressure (psi) 

Total dynamic axial load 

Local overpressure normal to the exterior face of a spherical 

dome (psi) 

Equivalent concentrated load for equivalent system (kips) 

Elastic dynamic buckling load (kips) 

Equivalent load on an element as a function of time 

Average value of P(t) for the far (or leeward) side of the arch 



l refl 



r-a 



'sb 



side- 



so 



soi 



stag 



Local overpressure on the front face of a buried 'rectangular 
structure (psi) 

Peak underground overpressure resulting at & given location from 
an underground burst (psi) 

Total vertical load (blast plus static) on column 
Force acting at any time, t 

Average value of P(t) for the near (or windward) side of the arch 
(psi) 

Dynamic plastic axial load capacity of column 

Reflected shock wave overpressure for angle of incidence of zero 
degrees (psi) 

Reflected shock wave overpressure for angle of incidence other 
than zero degrees (psi) 

Overpressure existing in the incident shock wave for any value of 
t-t d (psi) 

Overpressure existing in the incident shock wave when t-t , = t. 
(psi) d ^ 

Local overpressure on the exterior side walls of a rectangular 
structure 

Initial peak incident overpressure (psi) 

Peak overpressure of the shock wave -formed in the interior of a 
structure with openings (psi) 

Stagnation overpressure; the overpressure existing in a region in 
which the moving air has been brought completely to rest (psi) 

Axial load on the columns of a frame due to the vertical live and 
dead loads (kips) 

XIV 
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|r.P ) 
* V/ 



The time average of the total vertical load on the column in the 

time interval t < t < t 

ge gm 

P. . Average overpressure on the exterior back wall of a rectangular 
structure (psi) 

P. , Average net overpressure acting on back wall of a rectangular 
structure (psi) 

*, . J Peak value of the average overpressure on the exterior back 
' wall of a reci^ngular structure (psi) 

Average overpressure on the closed ends of a cylindrical segment 

/ . \ 

(psi) 

P. . Average net overpressure acting on front wall of a rectangular 
structure (psi) 

P~:fro . Average overpressure on the exterior of the front wall of a 



P~, - . Average overpressure on the interior front wall of a rectangular 
structure with openings (psi) 



rectangular structure (psi) 

. Average overpressure on t 
structure with openings ( 

P. fl Reflected shock wave average overpressure in the interior of a 
rectangular structure with openings (psi) 

P~i-roof Average overpressure on the interior roof of a rectangular struc- 
ture with openings (psi) 

P. . . Average overpressure on the interior side walls of a rectangular 
structure with openings (psi) 

P . Net average horizontal overpressure exerted on a rectangular 
structure (psi) 

P_ !. Average net overpressure acting on the roof of a rectangular 
structure (psi) 

P~roof Average overpressure on the front slope of a gable roof (psi) 

Average overpressure on the exterior of a rectangular structure 
(psi) 

P"roof Peak average overpressure on the front slope of the first gable 
of a multi-gabled roof (psi) 

? net Average net overpressure acting inward on sidewall of rectangular 
structure (psi) 

P . . Average overpressure on the exterior side walls of a rectangular 
structure (psi) 

P., Axial load determined by the intersection of the two linear 
portions of the P vs M curve of a steel beam section (kips) 

P(t) Actual load on a structural element as a function of time 

p /?atio of tensile reinforcement in reinforced concrete members to 
concrete area, A /bd 

Uniformly-distributed load intensity 

XV 
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TJ Critical steel ratio for reinforced concrete member 

^o 

p Equivalent static load for an arch (psi) 

p ' Ratio of volume of spiral reinforcement to the volume . of the 

concrete core (out-to-out of spirals) of a spirally reinforced 

concrete column 

Ratio of compressive reinforcement in beams to the concrete area, 

A^/bd 

Q Statical moment of the section about the centroidal axis 
Strength of a shear connector in composite construction 

Radiant energy on a unit area (cal/cm J 

q Unit drag pressure produced by the incident shock wave (psi) 
q. Dynamic design bearing pressure on soil I kips/ft y 
a Maximum bearing pressure on soil (kips/ft c ' 1 

q Maximum unit drag pressure produced by the incident shock wave 
(psi) 

R Radius of a spherical dome 

Total resistance of structural element or structural system (kips-) 
Dosage of gamma radiation (roentgens) 
Crater radius (ft) 

R. Total resistance of "A" portion of two-way slabs (kips) 
R Crater radius (ft) 

L> 

R C Horizontal static load resistance of shear wall at first cracking 
Horizontal dynamic load resistance of shear wall at first cracking 

Reynolds number of the high velocity wind in the incident shock 

wave 

Resistance of the equivalent system (kips) 

.:-implified notation of R the resistance function developed 

6 jo"- 1 - 

between the g and (g-l) floors of a multi-story building 
The average value of R in the time interval from to t' 

Resistance function developed between the g th and the i th floors 
of a multi-story building (kips) 

* Maximum resistance of the g th story of a multi-story building 
(kips) 

Resistance developed in the spring which connects m to the ground 
.^ Total resistance acting on the g th floor mass (kips) 
R m Maximum resistance developed by a structural system (kips) 

^ Plastic resistance of the "A" portion of two-way slabs where the 
edge is fixed (kips) 



dc 



R 



gi 
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R Plastic resistance of the "B" portion of two-way slabs vhere the 
^ edge is fixed (kips) 

R Maximum resistance of the equivalent system (kips) 
me 

R Fictitious maximum resistance (kips) 
mi 

R Resistance of a structural element or system at time, t (kips)' 

R Horizontal ultimate load resistance of a shear wall (kips) 

R Yield resistance of the structure (kips) 

R I Resistance in the elastic range (kips) 

R Maximum resistance in the elastic range (kips) 
1m 

R.. fl Maximum total resistance in the elastic range - "A" portion of 
two-way slabs (kips) 

R.. Maximum total resistance in the elastic range - "B ir portion of 
two-way slabs (kips) 

R(t) Resistance as a function, of tine 

R (t) Resistance of the g story column R as a function of time 

& 

(tj~j Tip** variation of the resistance of the equivalent single- 
g J e 1 degree-of -freedom dynamic system for the g tn story 

R Dosage of nuclear radiation without shielding (roentgens) 
R(x) Resistance as a function of displacement 

r Radius of gyration of section (in.) 

Resistance per unit length of a beam or per unit area of slab 
Ratio of web reinforcement * A /bs 

Distance from point of explosion (ft) 

Roentgens 

Ratio of steel reinforcing placed perpendicular to the steel beam 

(in composite construction) in excess of that required to carry 

the slab bending stresses 

R Average value of the resistance of the g story when the relative 
S displacement between the g tn and the g-l^h story is negative 

S Section modulus 

SE Strain energy absorption (kip -ft) 

(SE) Strain energy absorption of actual structure (kip-ft) 

& 

(SE) Strain energy absorption of equivalent structure (kip-ft) 
S 1 Section modulus about the weak axis 

s Coordinate axis at ^5 angle to x axis in x-y plane 

Spacing of stirrups and spacing of ties in reinforced concrete 
members (in.) 
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s Distance along s axis to centroid of loading (ft) 

s Distance along s axis to centroid of inertia force (ft) 

T Duration of the external load (sec) 
Resultant tension force (kips) 

T Fundamental period of vibration for complete spherical shell of 
c centroid thickness (sec) 

T Natural period of oscillation or fundamental period (sec) 



n 



Lovest natural period of circular arch with pinned or fixed ends 



.. 
T Natural period of the i mode of oscillation (sec) 

T Time or rise of a step-load (sec) 

T i; To---, T Time duration of the external loads on the first, second, 

o -*-"U 

and g floors of a multi-story building (sec) 

T 15C' T 2X"' T zX Time duration when tne effective loads on the first, 
IX dx. gx secOT x> Li an a gth stories are negative (sec) 

T-*- Period of oscillation of a fictitious system representing the 
^ g** 1 story (sec) 

t Thickness of concrete slabs (in.) 
Thickness of flange of beams (in.) 

Time measured after the arrival of the incident shbck wave (sec) 
Time variable 
Thickness of deep beam web 

t Time of arrival, or time required for the shock wave to travel 
from the point of explosion to the chosen location (sec) 

t Average flange thickness of standard steel beam ( in. ) 

t. Time required for overpressure on the rear face of a closed 

rectangular structure to rise from zero to its maximum value (sec) 

t fi Time required to clear the front face of a structure from the 
reflection effects (sec) 

t^ Time displacement factor; the time required for the shock front 
to travel from the frontmost element of a structure to the point 
or surface under consideration (sec) 

t g Time at which the limiting elastic deflection is reached 
t f Flange thickness of a WF steel beam (in.) 

t Time at which the limiting elastic deflection of the g th story 
is reached 

t m Time required for the vortex generated at the front face of a 
structure to travel a distance I, 1 across the structure (sec) 
Time required for maximum displacement of element or structure 
to occur (sec) 

t Q Duration of the positive phase of the incident shock wave (sec) 

XVIII 
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t Time of rise; time required for the overpressure in the incident 
shock wave to rise from zero to its maximum value (sec) 

t Thickness of stiff ener (in.) 

8 

t. Time required for the shock to pass over the length of the 
structure (sec) 

t Time of rise of an underground overpressure pulse (sec) 
t Web thickness of a steel beam or channel (in.) 
t Time required to reach yield point of material (sec) 

t.., tp...t , t .. Time sequence 

t., t_...t Time at which the maximum absolute displacements of the 
first, second and g** 1 floor masses are reached 

V Time required for the overpressure on the exterior rear face of 
a rectangular structure with openings to rise from zero to its 
maximum value (sec) 

t 1 Time required to clear the front face of a structure with openings 
from reflection effects (sec) 

t', t f Time at which the relative displacements X~ and X , respectively, 
are equal to zero 

t- Time immediately before t 
n n 

t+ Time Immediately after t 

Tiffle lag 
At Time interval used in numerical analysis 

At t , - t A time interval 
n n+1 n 

U Velocity of the incident shock front (fps) 

U. Velocity of the shock front of the shock wave formed in the 
interior of a rectangular structure with openings (fps) 

u Particle velocity in the incident shock wave (fps) 

Bond stress per unit of surface area of bar in reinforced con- 
crete design (psi) 

V Dynamic reaction (kips) 
Total shear (kips) 
Total vertical force on piles 

V. Total dynamic reaction along one edge "a," two-way slabs (kips) 

A 

V Average vertical shear in length, L 1 (kips) 

V_ Total dynamic reaction along one edge "b," two-way slabs (kips) 
B 

V Total column load in flat slab design (kips) 
c 

V Maximum shear capacity of web of deep beam section 
m 

V Maximum vertical shear force (kips) 
max 
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"dy 



m 



v n 



V 

v y 

v(t) 
w 



W a 
W 

V 

w 
wft) 

w 

v l 
X 

gm 



Vortex velocity (fps) 
Shear stress (psi) 

Dynamic shear yield strength 
: v (t \ Velocity of the g th floor at time, t n 

el n / 

Maximum shearing stress (psi) 
r(t \ Velocity at time, t n 

Velocity at t+ 
Velocity at t- 
Horizontal velocity of axis of rotation "0" 

Velocity of mass at time t g 

Static shear yield strength of steel (psi) 
Average unit shear stress in concrete (psi) 

Velocity as & function of time 

Total energy yield of an atomic bomb expressed KT of TNT 

required for an equivalent total energy yield 

Weight (ibs) 

Total load on element of structure (kips) 

Work done (ft-lbs) 

Dynamic peak blast load 

Work done on the actual system (ft-kips) 
Work done on the equivalent system (ft-kips) 
Maximum work done on the g floor mass (ft-kips) 
Absolute maximum work done on the g floor mass (ft-lbs) 

Maximum work done on the equivalent system by the equivalent load 
Fictitious maximum work done OD the equivalent system 

Work done as a function of time 

\ 

Uniformly distributed load (kips/ft) 
Length of channel shear connector ( in. ) 
Width of front face of rectangular building (ft) 

Weighting factor for vertical and batter piles 

Relative displacement in a story of a multi-story building (ft) 

Maximum relative displacement in the g th story (ft) 

Relative displacement of the mass g when the average external 
loads are applied statically 



I 



f g Relative displacement in the first, second, and g 
of a. multi-story building (ft), X = x - x 

XX 



th 



stories 
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Limiting elastic relative displacement for the g story column 



x Distance the incident shock wave travels after impinging on the 
frontmost element of a cylindrical segment or spherical dome to 
any point being considered 

Deflection of a structural element or system (ft) 
Distance of .any pile from the centroidal axis 

x Collapse deflection 

x Limiting elastic deflection (ft) 

x_ Forced solution of a dynamic system for a given external load 

x Absolute displacement of the g floor at time t 

x Maximum absolute displacement of the e floor (ft; 
gm jf o v 

x .., x , x .. Absolute displacements of the g floor at time, 
g,n-i g,n g,n+i 

n-l j n' ml 
x Maximum displacement (ft) 
X ,, x , x 1 Displacement at time t , , t , and t 1 , respectively 

x Initial displacement (ft) 

Horizontal displacement of axis of rotation "0" 

.c Displacement of an elastic system subjected to the peak load B 
acting statically (ft) 

xft ) Displacement of a mass at time t 

* ' e th 

x.., x 2 , x Absolute displacement of the first, second, and g floors 

8 of a multi-story building (ft) 

x Distance to centroid parallel to x axis (ft) 
X Q Horizontal velocity of axis of rotation "0" 

x" Horizontal acceleration of axis of rotation "0" 

.. / 2 i 

x Acceleration in the x direction ^ft/sec J 

y Displacement or deflection of a structural system (ft) 

y Deflection of actual element (ft) 

y Midspan deflection of actual element (ft) 

y- Midspan deflection (ft) 

y Limiting elastic deflection (ft) 

Deflection of equivalent system (ft) 

y Limiting deflection in the elasto-plastic range (ft) 

y Maximum displacement (ft) 

y Displacement at time, t (ft) 
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y 2m 

y 
y 

y 

y c i 

y 
z 
z 1 

z 
a 



a 



Vertical distance from axis of rotation "0" to mass centroid of 
rotating mass, m , of structure 

Deflection at midspan of column strip - flat slab (ft) 
Deflection at midspan of middle strip - flat slab (ft) 1 
Maximum elastic range deflection (ft) 
Maximum elasto-plastic range deflection (ft) 

Distance to the centroid parallel to y axis (ft) 

Vertical distance from axis of rotation "0" to mass centroid of 

total moving mass, m 

Acceleration in the y direction \ ft/sec j 

Velocity in the y direction (ft/sec) 

Velocity of displacement of the center of the beam (ft/sec) 

Increment of displacement (ft) 

Plastic modulus of the cross section 

Plastic modulus of the cross section in the weak direction 

Depth of the resultant compressive force in concrete tee beam 

Angle of incidence between the normal to the surface and the 

direction of propagation of the blast wave (degrees) 

Central angle of the arch 

Deflection coefficient for two-way slabs 

Design load ductility reduction factor 

Angular acceleration of structure about axis of rotation "0" 
(radians/sec 1 

Dimensionless ratio - 0.5P. 
Ductility ratio 

Maximum bending deflection due to the application of a unit load 
(ft/kip) 

Shear vail lateral deflection at first cracking 

Sub -area clearing factor for rectangular structures with openings 

Maximum shear deflection due to the application of a unit load 
(ft/kip) 

Shear wall lateral deflection at ultimate resistance 
Deflection at theoretical yield (ft) 

Angle of inclination to the horizontal of a gabled roof (degrees) 
Parameter used to define the location of a point on the surface 
cf .T cylindrical ce~nent or n spherical dome (degrees) 
End rotation of structural element 
Angular displacement 



. 
so 
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flL Joint rotation at "bottom of column 

D 

G C Rotation of beam at midapan (radians) 

Op Angular displacement of structure about axis of rotation "0" 

0_ Joint rotation at top of column 

/ \ 

V Viscosity of air in the incident shock wave ( Ib-sec/f t 2 ) 
Coefficient of friction 

V Poisson's ratio 

p Mass density of air in the incident shock wave ( Ib-sec /in. ) 

/ 2,4 2, 1A 
P B Mass density of soil ^Ib-sec /in. , or Ib-sec /ft J 

o cr Critical buckling stress (psi) 

Parameter used to. define the location of a point on the surface 
of a spherical dome 

Phase angle in the transient response (radians) 
Internal friction angle (degrees) 

0. Phase angle of the i mode (radians) 
7 Unit weight of soil (ibs/cu. ft) 

w Angular velocity of structure about axis of rotation "0" 
(radians/sec ) 
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Change 2 

ENGINEERING AND DESIGN 

DESIGN OF STRUCTUBBS TO BBSIST THE BPFBCTS OP ATCHIC WEAPONS 
PRINCIPLES OP DYNAMIC ANALYSIS AND DESIGN 

INTRODUCTION 

5-01 H1RPOSE AND SCOPE. This manual is one in a series issued for the 
guidance of engineers engaged in the design of permanent type military 
structures required to resist the effects of atomic weapons. It is appli- 
cable to all Corps of Engineers activities and installations responsible 
for the design of military construction. 

The material is based on the results of full scale atomic tests and 
analytical studies. The problem of designing structures to resist the ef- 
fects of atomic weapons is new and the methods of solution are still in the 
development stage. Continuing studies are in progress and supplemental 
Material vill be published as it is developed. 

The methods and procedures were developed through the collaboration 
of many consultants and specialists. Much of the basic analytical work was 
done by the engineering firm of Ammann and Whitney, New York City, under 
contract with the Chief of Engineers. The Massachusetts Institute of Tech- 
nology was responsible, under another contract with the Chief of Engineers, 
for the compilation of material and for the further study and development 
of design methods and procedures. 

It is requested that any errors.- and deficiencies noted and any sug- 
* gestions for improvement be transmitted^ HQDA (DAEN-MCE-D) WASH DC 20?U. 

5-02 KEKKKKNCES. Manuals - Corps of Engineers - Engineering and Design, 
containing interrelated subject matter are listed as follows : 

DESIGN OF STRUCTURES TO RESIST THE EFFECTS 
OF ATOMIC WEAPONS 

W 1110-3^5-1*13 Weapons Effects Data 

EM 1110-3^5-4 ill Strength of Materials and Structural Elements 

B* 1110-31*5-1*15 Principles of Dynamic Analysis and Design 
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KM 1110-31*5-1*16 
EM 1110-31*5-1*17 
EM 1110-31*5 -^IS 
EM 1110-31*5-1*19 
EM 1110-31*5-^20 
EM 1110-31*5-1*21 



5-OSa 



Structural Elements Subjected to Dynamic loads 

Single-Story Frame Buildings 

Mult i -Story Frame Buildings 

Shear Wall Structures 

Arches and Domes 

Underground Structures 



a. References to Material in Other Manuals of This Series. In the 
text of this manual references are made to paragraphs, figures, equations 
and tables in the other manuals of this series in accordance with the 
number designations as they appear in these manuals. The first part of the 
designation which precedes either a dash, or a decimal point, identifies a 
particular manual in the series as shown in the table following. 
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Items in the bibliography are referenced in the text by numbers inclosed in 
brackets. 

c. List of Symbols. Definitions of the symbols used throughout this 
manual series are given in a list; following the table of contents. 
5-03 RESCISSIONS. (Draft) EM 1110-31*5-1*15 (Part XXIII - The Design of 
Structures To Resist the Effects of Atomic Weapons, Chapter 5 - Principles 
of Dynamic Analysis.) 

5-<* BASIC PRINCIPLES USED DJ DYNAMIC ANALYSIS AMD DESIGN. Before discuss- 
ing the fundamental principles of dynamic analysis, it will be helpful to 
review briefly the principles used in the analysis of structures under static 
load. Two different methods are used either separately or concurrently in 
static analysis; one is based on the principle of equilibrium and the other 
is based on work done and internal energy considerations. 

Under the application of external loads, a given structure will be 
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deformed and internal forces vill be set up in its members . In order to 
satisfy static equilibrium, the vector sum of all the external and internal 
forces acting on any free body portion of the structure must be equal to 
zero. For the equilibrium of the structure as a whole, the vector sum of 
all the external forces and the reaction of the foundation must also be 
equal to zero. This principle is used in the analysis of statically de- 
terminate structures. 

The method based on work done and energy considerations is sometimes 
used when it is necessary to determine the deformation of a structure. In 
this method, use is made of the fact that the deformation of the structure 
causes the point of application of the external load to be displaced. The 
force then does work on the structure. Meanwhile, because of the structural 
deformations, potential energy is stored in the structure in the form of 
strain energy. By the principle of conservation of energy, the work done 
by the external force and the energy stored in the members must be equal. 
In static analysis, simplified methods such as the method of virtual work 
and the method of the unit load are derived from the general principle of 
energy conservation. 

In the analysis of statically indeterminate structures, in addition 
to satisfying the equations of equilibrium, it is necessary to include a 
calculation of the deformation of the structure in order to arrive at a 
complete solution of the internal forces in the structure. The methods 
based on energy considerations, such as the method of least work and the 
method based on Castigliano ' s theorems, are generally used. 

In the analysis of structures under dynamic loading, basically the 
same two methods are used, but the load changes rapidly with time and the 
acceleration, velocity, and hence, the inertial force and kinetic energy 
are of magnitudes requiring consideration. Thus, in addition to the in- 
ternal and external forces, the equation of equilibrium includes the in- 
ertial force, and the equation of dynamic equilibrium takes the form of 
Newton's equation of motion, which is given by equation (5-1)- 

(mass) (acceleration) = external force - internal force (5-1) 
As for the principle of conservation of energy, the work done must be equal 
to the sum of the kinetic energy and the strain energy, namely: 
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work done = kinetic energy + strain energy (5-2) 
The strain energy in equation (5-2) includes both the reversible elastic 
strain energy and the irreversible plastic strain energy. Thus the differ- 
ence between structures under static and dynamic loads is the presence of 
inertial force [(mass) (acceleration)] in the equation of dynamic equi- 
librium and of kinetic energy in the equation of energy conservation. Both 
terms are related to the mass of the structure. Hence, the mass of the 
structure becomes a very important consideration in dynamic analysis. 

In the above discussion, the two methods of analysis are described 
as if they were independent of each other. This is not the case. The force 
and energy equations are inter -dependent and are derived from each other. 
Although these two methods are not independent, the convenience of applying 
them to a given problem varies greatly depending on the particular problem 
under consideration. Generally speaking, the force equation is convenient 
for analysis of structures, and the energy equation is convenient for design. 
5-05 DYHAMICAmf EQUIVALENT SYSTEMS. In dynamic analysis, there are only 
three quantities to be considered. These are, (l) the work done, (2) the 
strain energy, and (3) the kinetic energy. To evaluate the work done, the 
displacement at any point on the structure under external distributed or 
concentrated load must be known. The strain energy is equal to the summa- 
tion of strain energies in all the structural elements, which may be in 
bending, compression, shear, or torsion. The kinetic energy involves the 
energy of translation and rotation of all the masses of the structure. The 
actual evaluation of these quantities for a given structure under dynamic 
load would be complicated. However, for practical problems this can be 
avoided by using appropriate assumptions. 

In order to simplify the problem, a giyen structure is replaced by 
a dynamically equivalent system. The distributed masses of the given 
structure are lumped together into a number of concentrated masses. The 
strain energy is assumed to be stored in a number of weightless springs 
which do not have to behave elastically. Similarly, the distributed load 
is replaced by a number of concentrated loads acting on the concentrated 
masses. Hence, the equivalent system consists merely of a number of con- 
centrated masses joined together by weightless springs and subjected to 
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concentrated loads which vary with time. This concentrated mass -spring- 
load system is defined as an equivalent dynamic system. 

The reduction of a given structure to an equivalent dynamic system 
involves the principle of dynamic similarity which is merely the require- 
ment that the work done, strain energy, and kinetic energy of the equiva- 
lent system must be identical respectively with those of the given struc- 
ture [2, 7] The details of this principle and illustrative examples are 
given in paragraphs 6.13 through 6.2k in EM 1110-3^5-^16. 

For the time being it is sufficient to state that any given struc- 
ture, be it a multi-story building or a structural member such as a beam or 
column, can be approximated by an equivalent dynamic system for the purpose 
of dynamic analysis. The discussion in this manual is limited to the 
analysis and design of equivalent dynamic systems consisting of concentrated 
masses, springs, and loads. 

a. Basic Dynamic System. The simplest dynamic system consists of a 
concentrated mass supported by a weightless spring and subjected to a con- 
centrated load. This is shown in figure 5-1- 



f(t) 




R 





Figure 5.1, Basic dynamic system 



In figure 5-1, the following notation is used: 
f(t) - external applied load which varies with time 
m mass 
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R m internal resisting force which depends on the displacement 

x displacement which also varies with time. 

The directions in which the variables f(t), R and x are considered to 
be positive are indicated by arrows in figure 5.1. 

Suppose at a given time, t , the displacement, velocity, and ac- 

p O 

celeration of the mass are x , dx/dt, and d x/dt respectively; then, 
the following expressions may be written: 



Internal resisting force = R * kx when the spring is elastic 



Inertial force = m 



d 2 x 



dt 
x t 



Work done by the load f f(t) dx = f f(t) f at. 

J 

v2 
Kinetic energy " f m ( '3T , 

/x 12 

Rdx = ^ kx when the spring is elastic 
^-L 



Bewton's equation of motion is: 

d 2 x 
d? 

and the equation of conservation of energy is 



(5-3) 



Rdx 



These are the two basic equations used in dynamic analysis. Before attempt- 
ing to evaluate the displacement, x , from these equations for' a given ex- 
ternal load, f(t) , a few concepts about the external load and the resist- 
ance function need clarification. 

(1) External Load and Work Done. load Duration and Load Impulse. 
A typical load applied to a dynamic system is shown in figure 5.2. This 
load acts over an interval of time, T , which is defined as the duration 
of the load. The impulse of the external load, denoted by H , is equal 

6 
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to the area under the load-time curve; namely: 

T 
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H 



f(t)dt 



(5-5) 



Expression for Work Done. The work o 
done by the external load up to any time, t , 
is given by: 




T 
It (56) Time 

Figure z.2. A. typical 
dynamic load 

In order to evaluate the work done, it is neces- 
sary first to determine the velocity, dx/dt- , which can be obtained from 
equation (5-3) by integration. This is given in symbolic form as: 



[f(t)-]dt 



(5.7) 



Substituting equation (5-7) into equation (5-6), the expression for work 
done then is, 



dt 



(5.8) 



Equation (5-8) indicates that the work done by the load varies with time. 
A typical case is shown in figure 5-3&- At time, t , when the displace- 
ment of the mass reaches its maximum value, the work done is also a maximum 

The maximum work done W by the external load is obtained when the in- 

m 

tegral limits in equation (5-8) are from to 

,-t 



m 



W 



m 



[f(t) - R] dt 



dt 



(5-9] 



Because of the presence of m and R in the integrand, it is evident that 
the ^x^"" work done depends not only on the external load but also on the 



mass and resistance function, R , of the dynamic system. 

Equation (5.9) is the general expression for the maximum work done. 
For the cases shown in figures 5-3b and 5-3c, the external load is zero 
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over the time interval from T to 

t . The value of W is a constant 
m 

after T . Hence, for the case where 

t > T , the maximum work done is 

m ' 

given by: 



f m T 



[f(t) - R] dt >dt (5-10) 

Pure Pulse Load. For the 
case shown in figure 5-3c, the value 

of t is much greater than T and 
m 

the value of internal resisting 
force, R in equation (5.10) is 
negligible compared with the value 
of the external load over the entire 
load duration, T . Under this con- 
dition, the external load is defined 
as a pure pulse load. In equation 

(5-10), if it is assumed that the resistance R is zero for the full load 
duration, the maximum work done by a pure pulse load is given by: 

ft 



Figure 5.3. External load, displacement function 
and work done curves 



W = / 

m J 



f(t) 



f(t) dt 



dt * 



(5-11) 



where 



H I f(t)dt = the impulse of the external load. 
'O 



A pure pulse load corresponds to an impact load. During impact, the 
internal resisting force and the strain energy are both assumed to be zero. 
After impact, the mass, m , acquires an initial velocity equal to H/m . 
The kinetic energy of the mass after impact is given by equation (5.11) 
which is also the maximum work done by the load. This equation indicates 
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that for a pure pulse load, the maximum work done depends on the load Im- 
pulse, but is independent of the detailed shape of the external load. 

Fictitious Maximum Work Done. In many actual problems, the in- 
ternal resisting force, R , in the time interval from to T is not 
small enough to be. neglected. From equations (5-3) and (5.10), it can be 
seen that if the internal force in the time interval from to T is 
neglected, the actual acceleration, velocity, and displacement will always 
be larger than the corresponding values when the resistance is considered. 
Likewise, the internal force in equation (5-9) always acts to reduce the 
work done by the external load. In other words, the actual maximum work 
done on a dynamic system will always be smaller than the maximum work done 
on the same system wher.e the load is considered to be a pure pulse load. 
Hence, the work done as computed from equation (5-11) is defined as the 
fictitious maximum work done, namely, 



The maximum work done, W , as computed from equation (5-9) is al- 
ways smaller than W . The ratio W /W depends on tn , R , and on the 
detailed shape of the external load-time curve. The work done ratio is an 
import"- 4 - non-dimensional parameter to be considered in the design of dy- 
namic systems. 

(2) Spring Constant, Resistance and Strain Energy. From para- 
graph 5-05a, it is seen that the function of the internal force is to re- 
sist the movement of the mass, hence the internal force is defined as the 
resistance of the structure. The variation of the resistance with dis- 
placement of the mass is expressed by a resistance function. 

In figure 5.1, the strain energy is assumed to be stored in a spring. 
For a linearly elastic spring, the resistance function is equal to kx , 
and, the strain energy is ~ kx 2 . For the general case, the strain energy 
is equal to jC X Rdx . Then the spring constant, resistance function, and 
strain energy are related to one another by the following two equations: 

. R (5-13) 
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(5.110 



where 

SE - internal strain energies which include both the reversihle 

elastic strain energy and the irreversible plastic strain energy 

R resistance 

k * spring constant 

Equations (5-13) and (5- 1*0 * obtained from the basic dynamic system shown 
in figure 5.1. For complicated cases when more than one displacement is 
necessary to specify the configuration of the dynamic system, these two 
equations still hold true except partial derivatives are used. 

The relationship between resistance and strain energy as expressed by 
equation (5-13) indicates that two different methods may.be used to evaluate 
the resistance function. One method is based on the direct evaluation of 
internal resisting force. The other method is based on the evaluation of 
strain energy, and the resistance is then derived from equation (5.13). 

The resistance when plotted against the displacement, x , usually 
takes one of the forms shown in figure 5-^a. For structures made of brittle 
mterials, the resistance function is generally as indicated by curve A. 
For structures made of ductile materials with marked yielding, the 
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Figure 5.4. Resistance functions 
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resistance function is given by curve B. For curve C, the slope of the re- 
sistance function is negative for large values of x . This is generally 
not a property of the material used, but is due to: (1) the combination of 
axial and lateral load on a given member of a given structure, and/or (2) 
the effect of elastic or inelastic instability. In evaluating the response 
of structures under dynamic load, if better accuracy in the result is de- 
sired, the resistance function as shown in figure 5. he. should be used. But 
the method of preliminary design is usually based on a number of design 
charts which are evaluated on the basis of idealized resistance functions 
such as shown in figure 5-^t> where one or two straight lines are used to 
approximate the actual resistance function. The type of idealized resist- 
ance to be used is discussed in paragraph 5-07- 

b. Degree of Freedom of Dynamic Systems. The basic dynamic system 
considered in the preceding paragraph consists of one mass, one spring and 
a concentrated load. Hence, a single displacement variable, x , is suf- 
ficient to describe its motion. This basic system is called a single- 
degree-of -freedom dynamic system, or a single-degree dynamic system. 

For the case of a multi-story building, more than one displacement 
variable is needed to describe its motion. The equivalent system consists 
of several' concentrated masses connected by springs. 

The degree of freedom of a dy- 
namic system is defined as the num- column, 
ber of independent displacement vari- 
ables needed to specify completely 
the configuration of the system. 
Based on this definition, the dy- 
namic systems shown in figures 5-5& 
and b have single, and two degrees 
of freedom, respectively. 



Eiternol 
Lood 



, 





c 



For ordinary building struc- 
tures under lateral loads, the mass 
of the structure can be assumed to 
be concentrated at the floor levels 
and at the roof. The framework be- 
tween floors can be assumed to be 
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Figure 5.5. Single- and two- degree- of- freedom 
dynamic systems 
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weightless splines or flat springs. The degree of freedom of a building is 
therefore equal to the number of stories of the building. In certain spe- 
cial cases, the displacements of any two floors can be assumed to always 
bear a constant ratio. The degree of freedom is then reduced. A multi- 
story building may even be approximated by a single-degree dynamic system. 
However, the error involved in this simplification may be very large. It 
is not recommended unless it is only for a first approximation in analysis 
or design. 

Beams, columns, slabs, etc., which have a continuous distribution of 
mass as shown in figure 5-5& are actually infinite-degree dynamic systems. 
However, if it is assumed that the deflected shapes of the member at any 
two instants of time are geometrically similar to each other, only one dis- 
placement is necessary (for example, the deflection at midspan). Thus, an 
inf inite-degree-of -freedom system is reduced to a single-degree dynamic 
system. The method of evaluating the parameters of the equivalent dynamic 
system for these members is discussed in paragraphs 6-13 through 6-2U, 
EM 1110-3^5-Ul6 . 

c- Difficulties in Application of Dynamic Systems. A basic dif- 
ficulty in dynamic analysis results from the fact that the method of super- 
position is not directly applicable. If there are several loads applied 
simultaneously on a dynamic system with elastic springs, the maximum values 
of the displacements produced by individual loads do not generally occur at 
the same time. The maximum value of the displacement produced by all the 
loads is not equal to the summation of the maximum values produced by in- 
dividual loads. Hence, in order to apply the method of superposition, the 
so-called phase differences in the dynamic responses to individual loads 
must be considered. Moreover, in some cases, the equivalent dynamic system 
of a structure may be non- linear, for example when a certain member is 
stressed beyond the yield point or when there is a vertical load as well as 
a lateral load acting on the structure. For non-linear dynamic systems, 
the method of superposition is not applicable at all. Hence, in general, 
all the loads acting on a dynamic system have to be considered 
simultaneously . 

In dynamic analysis, a given structure is approximated by its 
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equivalent dynamic system. Any structure is made up of a number of ele- 
ments such as beams, columns, girders, and structural connections. Each 
member is a distributed mass and possesses an infinite degree of freedom. 
Rigorously speaking, any structure is a multi-infinite degree dynamic 
system. However, by making certain assumptions and by applying the prin- 
ciples of dynamic similarity, a complicated structure is simplified to an 
equivalent dynamic system with a finite degree of freedom. Essentially, a 
distributed mass-load-spring system is replaced by a lumped mass-load- 
spring system. The accuracy of the results obtained from the equivalent 
system as compared with those obtained from the actual structure depends 
entirely on the accuracy of the assumption used in the derivation of the 
equivalent system. In some cases, it may involve large error. It should 
be understood that the replacement of an actual structure by an equivalent 
dynamic system is one of the basic approximations used in dynamic analysis. 

Even for those cases when the equivalent system is acceptable as a 
good approximation for the given structure, there are difficulties in the 
evaluation of its parameters. One typical example is the evaluation of the 
resistance of a one-story frame building. The internal force which resists 
the lateral movement of the roof is equal to the total shear force at the 
top of the columns and this shear force is developed because of the bending 
and shear deformations of the columns. The bending deformation, hence the 
resistance, depends on the degree of fixity of the ends as well as on the 
size of the columns. The determination of the end fixity involves the size 
of girders, foundations, and the type of end connections. In practical 
problems, for the sake of simplicity, the two extreme conditions of fixity 
are generally used, that is, the end is assumed either hinged or fixed. 
For the determination of the resistance function of a member with stresses 
exceeding the elastic limit, it is necessary to consider the variation of 
the elastic limit with the rate of strain which is discussed in paragraph 
4-03a, EM 1110-3^5-^13. There are still other effects to be considered 
such as the effect of dead weight and vertical load on the roof. Because 
of these factors, the actual resistance function is a complicated function 
of the displacement. Although in the analysis of structures, these factors 
may be partially taken into account in order to obtain a resistance function 
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which is a closer approximation to the actual resistance function, idealized 
resistance functions are always used in design. 

The same difficulty is encountered in the determination of the load as 
in the determination of the resistance function. The load is determined 
either directly or indirectly from the blast pressure which is generally a 
complicated function of time. Here again, the actual variation of load can 
be taken into account in analysis, but simplified loads are required for 

design. 

From the above discussion, it is seen that the parameters of the dy- 
namic system are subject to uncertainty in-their evaluation. The analysis 
and design of a structure based on its equivalent dynamic system can only 
be carried out approximately. This is especially true in design since 
idealized resistances and simplified loads are used. Although rules are 
given at appropriate places which will help the designer to make reasonable 
assumptions and approximations, a basic understanding of the effect of dy- 
namic loads is necessary for the intelligent application of these rules to 
actual problems. 

In paragraphs 5-06 through 5-08, methods of analysis of single- 
degree dynamic systems are discussed. The results of a systematic anal- 
ysis of single-degree dynamic systems for several simplified loadings are 
plotted in non-dimensional design charts in paragraph 5-10. The applica- 
tion of these charts to the design of a. single-degree dynamic system is 
discussed in paragraphs 5-12 through 5-1^ The analysis and design of 
multi-degree dynamic systems are discussed in the remaining portion of 
the manual. 

ANALYSIS OF SIWGLE-DEGREE-OF-FKEEDOM SYSTEMS 

5-06 INTRODUCTION. The analysis of single-degree dynamic systems consists 
of the evaluation of the displacement using Newton's equation of motion 
(equation 5.3). Two methods are generally used in solving this type of 
differential equation. One is the rigorous method in which the solution 
is obtained directly from the differential equation. This method is prac- 
tical only when both the load and the resistance function can be expressed 
in simple mathematical forms. The other method is the numerical method in 
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which the solution is obtained by a step-by-step integration. The numerical 
netfaod is generally applicable for any type of load and resistance function. 
The two methods of analysis are discussed in paragraphs 5-07 and 5-08 re- 
spectively. A numerical example is given in paragraph 5-09 to illustrate 
the different methods of analysis. 

5-07 ANALTSIS BY RIGOROUS METHOD, a. General. The problem considered in 
this paragraph is the evaluation of the displacement using Newton's equation 
of notion for a single-degree dynamic system as shown in figure 5.1. The 
equation of motion is: 

2 

mi-|.f(t)-R (5.15) 

d1T 

For a given load and a given resistance function, the procedure for 
evaluating the displacement is the same as that for solving a second order 
differential equation. First, the complementary function of the differen- 
tial equation is obtained by assuming the load to be zero. There are two 
constants of integration in the complementary function. Next the partic- 
ular integral for the 'given load is evaluated. The complete solution for 
the displacement is the sum of the complementary function and particular 
Integral which are defined respectively as the transient and the forced 
solutions in dynamic problems [2, 3, 7J- From the known initial condi- 
tions, that is, the displacement and the velocity of the mass at the 
initiation of load, the constants of integration are evaluated. Thus an 
expression of the displacement as a function of time is obtained. It is 
obvious that the displacement function depends on the mass, m , the dy- 
namic load, f(t) , and the resistance function, R . 

The rigorous method of analysis is generally impractical for problems 
when either the dynamic load or the resistance 'is a complicated function. 
However, many important concepts related to the basic principle of dynamic 
analysis and design can be illustrated when this method is applied to 
Idealized cases in which both the load and the resistance are of simple 
mathematical form. 

Three types of dynamic load of simple form are shown in figure 5-6 
together with the expressions for them. Each type of load is completely 
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(a) Rectangular 
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Figure 5.6. Three idealized dynamic loads 



Three types of idealized resistance functions as shown in figure 5-7 
are chosen for study. The dynamic systems characterized by these resistance 
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Figure 5.7. Three types of resistance function 

functions are defined as: (a) linearly elastic, (b) completely plastic, 
and (c) elasto-plastic . The expressions of these resistance functions are 
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also given in figure 5.7. It is noted that the linearly elastic and. com- 
pletely plastic systems are only special cases of the elasto-plastic system. 
The rigorous method of analysis of these basic types of single-degree dy- 
namic systems are given separately in the following three paragraphs. 

t>. Linearly Elastic System. The resistance f unction of a linearly 
elastic system is always equal to kx . Substituting kx for R in 
equation (5-15), the equation of motion becomes: 

- kx (5-16) 



at 



Referring to a textbook on Differential Equations, the complete solution of 
equation (5-l6) is given by: 

x . C sin2Tr jL - j + x (5-17) 



/2Tr jL 



where 

T = 2itVm/k = natural period of oscillation 

x = forced solution 
C = constant of integration 
= phase angle (constant of integration) 

The sinusoidal term of equation (5-17) is the transient solution which is 
associated with the oscillatory nature of the dynamic system. The two con- 
stants of integration, .C and , are determined from the requirement 
that the displacement and velocity of the mass must satisfy certain pre- 
scribed initial conditions. For structures under dynamic load, the usual 
initial conditions are: 

at t = 0, x = 0, and dx/dt = (5-18) 

Once C and are determined, the displacement at any time is obtained. 
Equation (5-17) shows that, for any given dynamic load, the speed of 
response of the linearly elastic system depends on the natural period of 
oscillation, T , in the transient solution. On the other hand, the 
rapidity of variation of the applied load depends on the load duration, T 
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Thus the non-dimensional ratio, T/T n , is expected to be an important 
parameter in dynamic analysis and design. 

AB an illustrative example, the case of a triangular external load 
given in figure 5.6b is considered. The expressions for the load are dif- 
ferent depending on whether t is greater or smaller than T . The cor- 
responding displacement functions are also different. 

For t < T , the expression for the load is: 

f(t) - B (l -|) (5-19) 

The forced solution for this load is: 

(5-20) 



, 2 (i - | 



Substituting equation (5-20) into equation (5-1?) and taking the initial 
conditions given by equation (5-18) into account, the arbitrary constants 
C and $ can be evaluated. These are given by: 




(5-22) 

n/ 

Thus, in the time interval from to T , the displacement is given by: 




x 

where 

n 

x s * k = sta ' tic deflection produced by the peak load, B . 
After t T , the external load is zero. The forced solution is also zero. 
Prom the displacement and velocity at time, T , determined from equation 
(5.23), the constants of integration of the transient solution after t > T 
can be determined. The displacement function after t > T is given by: 
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Equations (5-23) and (5-24) show that the non-dimensional displace- 
ment variable x/x at a given time, depends only on T/T . In figure 
5-8a, the displacement of a linearly elastic system subjected to a tri- 
angular load is plotted versus time 
for three values of T/T , namely, 
T/T 0.2, 0.5, and 2.0. The cor- 
responding cases of rectangular load 
and step load with rise time, T , 
are plotted in figures 5.8b and c. 
These curves show: (l) the maximum 
displacement, (2) the time when the 
maximum displacement occurs, and 
(3) that the time for one complete 
oscillation is equal to T Q , and 
(k) the displacement at any time. 

c. Completely Plastic System. 
The resistance of a completely 
plastic system is a constant equal 
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than the maximum elastic deflection 
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Figure 5.8. The displacement functions of elastic 
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and the system treated as a completely plastic system. 

Substituting R_ for R in equation (5-15), the equation of motion 



ffl 



of a completely plastic system is: 



(5-25) 



dt 



Once f(t) is given, the displacement function can be evaluated directly 

by integration. 

For the case of a rectangular load given in figure 5.6a, the expres- 

sions of the load are: 

f(t) = B f or t < T (5-26) 

f(t) = f or t > T (5-2?) 

Substituting equations (5-26) and (5-2?) into equation (5-25), integrating 

and taking equation (5-l8) into account, the displacement functions of the 

mass for a rectangular load are given by: 



for t < T 



, f . (l - ^)(i 



(5-28) 



f or T < t < t 



x 
F 



(5.29) 



where 



. 1 BT 

V ' 

x s " 2 



m 



and 



t = the time when the maximum displacement x is reached. 
m m 

Substituting equation (5-28) into equation (5-6) and integrating, 
the work done by the external load up to any time, t , which is less 
than T can be obtained. This is given by equation (5-30). 

.2 



where 



(5-30) 



o ~" ^ see et l ua ' t:l - on 5-12) 



1 H 2 
* o 

1 (BT) 2 
= 2 m 
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The maximum work done by the load occurs at time, T , and is given by: 

(5-31) 



i-ifl ^ 

\ n 1 B; 



The preceding four equations show that the non-dimensional quantity 
R /B is an important parameter in the analysis and design of completely 
plastic systems. In figure 5-9a, the displacement and the work done fcy the 

load are plotted versus time for a rectangular load with B = 2R . The 

m 

corresponding functions for a triangular load are plotted in figure 5-9b- 
These curves show: (l) the displacement function for a completely plastic 
system is non-oscillatory; (2) the maximum displacement; (3) the maximum 
work done; and (^) the time t when the maximum displacement x occurs. 




1.0 

Time Rotio 



1.0 
Time Ratio 



(a) Rectangular- Load (b) Triangular Load 

Figure 5.9. The displacement and work done on a completely plastic system B 2# m 

d. Elasto -Plastic System. This is the dynamic system which is com- 
monly used to represent a structure as a whole, as well as its members, for 
the purpose of analysis and design of blast resistant structures. Under 
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entirely in the elastic range. Many ductile construction materials are 
characterized by a large yield range which enables them to absorb a large 
amount of energy. Use of this property in design results in much lighter 
nentoers than would be required should the energy absorption ability of the 
yielding region not be utilized. 

For an elasto-plastic system, the resistance function is given by: 



R = kx for x < x. 



R . R. 

m 



for x > x 



(5-32a) 
(5-32b) 



After the maximum displacement is reached, the dynamic system rebounds and 
the resistance function is given by: 

X (5-32c) 



M*- 

R . - R 



f or 



for 



- ) 



( x m ' x ) 



> 2x 



(5-32d) 



m \ m / e 
For the rectangular and triangular loads shown in figure 5-6, the expres- 
sions of the load are different depending on whether t is greater or 
smaller than T . Hence, the six differential equations given below are 
needed to express the equation of motion. 



t < T, x < x 



t < T, x > x , m 




f(t) - kx 



f(t) - 



m 



(5-33a) 
(5-33t>) 
(5-33c) 
(5.33d) 
(5-33e) 
(5-33f) 



The first four equations express the motion before the maximum displacement 
is reached and the last two equations express the motion when the dynamic 
system is in elastic rebound. 

The sequence of using these equations for the evaluation of the 
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displacement function of a given problem depends on: (l) the time, t , 
when the maximum elastic deflection is reached, and (2) the time, t , 
when the maximum displacement is reached. 

For the linearly elastic system, T/T is an important parameter 
while R/B is important for a plastic system. Both parameters are im- 
portant for an elasto-plastic system. As an example, the case of a tri- 
angular load as shown in figure 5-6b is considered. The numerical values 
of the non-dimensional parameters are assumed for illustrative purposes as: 



c f = ~ 

n 



Equation (5-33a) is first used. Following the method described in para- 
graph 5-07b, the displacement function is given by: 



where 



B 



-1 T 
tan i 2i- 

n 




(5.3>0 



Equation (5- 3*0 is true so long as t < T and x < x . Substituting 



B 
R k 



into equation (5-3^), the time, t , when the maximum elastic deflection, 

C 

x , is reached can be evaluated. This is given by: 

t /T = 0.64 

The displacement and velocity of the mass at time, t g , are respectively 
given by: 



w 



x = C_ 
e R k 



a 1.11 



_ 

kT 



23 



EM 1110-31*5-1*15 
15 Mar 57 



5-07d 



In the time Interval, t < t < T , equation (5-33b) should be used. 
The displacement function is obtained by direct integration. This is given 
by: 



B 



C 

R 



RT 



2 r 



(5.35) 



At time, t T , the displacement and velocity are respectively given by: 

x(T) = 0.831 | 

v(T) = 0.702 (r 

In the time interval, T < t < t ffi , equation (5 -33d) should be used. The 
displacement ' in this time interval is: 



B 
x "k 



(5.36) 



Prom equation (5-3^)> the maximum displacement, x , is given by: 



x = 0.971 ^ 
m ^ k 



and the maximum displacement occurs at 



m 



-nF = 1.397 



After t > t m , the dynamic system is in elastic rebound, and equa- 
tion (5-33f ) should be used. The displacement after t = t is given by: 



*-! 



m 



0.1*71 -t- i cos 



(5-37) 



The displacement function as expressed by equations (5.3*0, (5.35), (5.36), 
or (5-37) is plotted in figure 5-10 in non-dimensional form. 
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Figure 5.10. Displacement and work done by a triangular loud an at tlaato-plastic system 



.---0.3, 



m 



Once the displacement function is known, the work done by the ex- 
ternal load at any time, as expressed by equation (56), can be evaluated 
and this is also plotted in figure 5-10- The maximum work done by the ex- 
ternal load, which is equal to the total strain energy at the maximum dis- 
placement, is given by equation (5-38). 

S 

(5-38) 

As shown in figure 5.10, the ratio of W m /W for the given problem is 
equal to O.SlU. 

5-08 ANALYSIS BY NUMERICAL METHODS, a. Need for numerical Method. The 
rigorous method described in paragraph 5-07 is practical only for the cases 
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f(t) 




5-08a 

when both the load and the resist-- 
ance are simple mathematical expres- 
sions. However, for an actual prob- 
lem, the load may be quite compli- 
cated. A typical case is shown in 
figure 5.11. Four separate equation 
are needed to express f(t) . For 
this load applied to an elasto- 
t plastic system, the equations of no- 
tion are expressed by eight separate 
differential equations. The solutio 
of these equations by rigorous metho 
is generally impractical. 

The rigorous method is further complicated if the resistance functior 
is not a simple shape as shown in figure 5-7- Two practical cases are 
shovn in figure 5-12. One is the idealized resistance function of a fixed- 
end beam, and the other is of the columns of a single story frame building 
when the vertical load is taken into consideration in the evaluation of . 
the resistance. 

For practical applications, when the load and resistance functions 
are complex, because of the impracticability of solving the rigorous 

R R 



Time 

Figure 5.1 1. A typical dynamic load 





Displacement 



Displacement 



(a) Resistance Function 
of a Fixed-end Beam 



(b) Resistance Function 
of a One-Story Building 

Figure 5.12. TKO practical resistance functions 
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differential equations, it is necessary to have other methods of analysis. 
Many methods have been developed. One is the method of numerical integra- 
tion which is discussed in this section. 

Perhaps it should be mentioned that these methods are by no means in- 
exact. Any degree of accuracy desired in the end result can be obtained. 
But, for the convenience of computation, these methods always embody sim- 
plifying assumptions vhich affect the accuracy of the results. The inac- 
curacies are not due to inherent properties of the numerical methods, but 
rather due to the simplifying assumptions used in carrying them out. 

b. Basic Principles of Numerical Analysis. The differential equa- 
tion of motion (equation 5-15) expresses implicitly the acceleration of a 
dynamic system at any time, a m " It is known that the velocity 
and displacement can be obtained 
from the acceleration by inte- 
gration. The numerical method 
of evaluating the displacement 
from equation (515) is called 
the method of numerical inte- 
gration. 

The load, resistance, ac- 
celeration, velocity, and dis- 
placement are plotted versus 
time in figure 5-13 for a typ- 
ical case. 

Suppose t , tp 

*2'"*n - 1' *n' *n +l'" is 
a time sequence. The time in- 
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f(t)-R 
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terval from t 
denoted by At 



to t. 



is 



The dynamic 



load is assumed to be initiated 



at t t 



The acceleration, 



velocity, and displacement at 
t fl are denoted by a Q , V Q , 
and x , respectively. If the 
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Figure 5.13. Load, resistance, acceleration, velocity 
and displacement versus time 
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acceleration in the time interval, At n , is represented by a(t ) , the 
velocity and displacement at t n + 1 are given respectively by the follow- 
ing equations: 



v n 



1 a(t) dt (5-39) 



These two equations indicate that the velocity and displacement at V ^ 
can be obtained by extrapolation from the corresponding values at t 
once the acceleration in the time interval, At , is known. 

In the analysis of structures under dynamic load, the velocity and 
displacement at the initiation of load are given as zero, that is: 

v - and x ( 

o o 

In applying equations (5-39) and (5-^0), the values of v., and x, can be 
obtained provided that a(t) is known in the time interval from t to 
t., . Prom the values of v.. and x, , the values of v_ and x~ can be 
obtained. This process can be continued until the values of v and x 
for any values of n are obtained. Hence, in the above numerical proce- 
dure the displacement, x, , x ? , x_ . . . , x . . . are evaluated by a step- 
by-step extrapolation method starting from the known initial conditions. 
Here the difference between the rigorous method and the numerical 
method can be perceived. In the rigorous method, the value of x is ob- 
tained as a continuous function of time while in the numerical method, the 
value of x is obtained at discrete values of time. In other words, the 
actual relationship between x and t is given by a differential equa- 
tion, but in the numerical method, this relationship is replaced by a 
finite difference equation. As the time interval in the sequence of time 
is reduced, the finite difference equation approaches the actual differ- 
ential equation as a limit. It can, therefore, be concluded that the ac- 
curacy of the numerical method based on the finite difference equation for 
most problems depends almost entirely on the magnitude of the successive 
time intervals. 
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Many extrapolation formulas have been derived for the solution of the 
finite difference equation [5] . For engineering application, simplified 
extrapolation formulas are often used which are obtained by assuming a 
simple acceleration-time relationship in any time interval, At Q . Two 
procedures are developed in paragraphs 5-08c and 5-08d. In one case the 
acceleration is assumed to be varying linearly with time. In the other 
case, the acceleration is assumed to consist of a series of pure pulses. 

c. Linear Acceleration Extrapolation Method. In this method, the 
acceleration, a(t) , in any time Interval, At , as shown in figure 5- 
la assumed to vary linearly with time as shown in figure ' 
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(b) Approximate Acceleration 
Curve 



Figure 5.14. Comparison of actual and approximate acceleration curve 
used in linear acceleration extrapolation method 

In the time interval, At n , a(t) is then given by: 



a(t) 



^e) ( - ) 



At 



n 



Substituting equation (5.1*2) into equations (5-39) and (5^0), the recur- 
rence formulas for velocity and displacement are given respectively by: 
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n 

n + 1 




n 4- 




Once a and a . are known, these tvo equations are applicable to the 
evaluation of the velocity and displacement at t .. from the known values 



at t 



n 



According to definition, a(t) is given by: 



dt 



* l) - 4" + 1 



are 



In equation (5-^5), fit ,) is known, "but the resistance, Rlt , J , 

\ n T J-/ \ n "r -L/ 

may depend on x , , which is still an unknown. Hence, in the applica- 
tion of equations (5-^3) and (5.kk), a trial and error procedure has to be 

adopted. A value of x T is assumed as a first trial, x , can 

n + 1 L / ,'n + l 

be assumed to be X Q + V Q [At n | from which Rlt J and a - 
obtained. The value of x , is computed from equation (5.W-) and the 
computed value is compared with the assumed value. This procedure is re- 

peated until the assumed and computed values of x ., agree with each 

n *t* x 

other. 

Of course, in the plastic region, the resistance is a constant; 
x , can be obtained directly from equation (5.UU) and the trial and 
error procedure is not necessary. 

There are two basic characteristics of this method. First, it is, 
In general, a trial and error method, hence it may be inconvenient to carry 
out, and secondly, the acceleration is assumed to vary linearly with 'time, 
thus the variation of acceleration in the time interval, At , is par- 
tially taken into account. So this method is called a trial and error 
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linear acceleration extrapolation method. 

In the application of the above method, the successive time intervals, 
At , At , At 2 , . . .At do not have to be equal although they are chosen to 
te equal for the convenience of computation. To determine the magnitude of 
the time interval to be used for a given problem, a number of factors have 
to be considered. First, the time interval must be chosen such that within 
It the acceleration can be reasonably approximated by a straight line. 
Secondly, the convergence of the trial and error procedure must be rapid. 
Dilrdly, the number of steps necessary for the evaluation of the maximum 
displacement must be reasonable. For the first and second considerations, 
the tine Interval should be small while for the third consideration, the 
tine Interval should be as large as tolerable. Because of these conflict- 
lag requirements, it is found that a time interval equal to about T _/6 is 
a good compromise to be used in the elastic region, but a larger time in- 
terval can be used in the plastic region provided that the external load in 
the tine Interval is a reasonably straight line. 

d. Acceleration Impulse Extrapolation Method. In this method, the 
actual acceleration curve shown in figure 5-15a is replaced by a train of 
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Figure 5.15. Acceleration impulse extrapolation method 
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t.., 



equally spaced Impulses occurring at t , t.. 
of the acceleration impulse at t is given by: 



t_, ...t 



5-08* 

The magnitude 



(At) 



(5.1*6) 



where 



At 



" ' * " 



This is shown in figure 5-15b. Since an impulse is applied at t n , there 
is a discontinuity in the value of velocity at t . In the time interval 

from t to t , the velocity is constant and the displacement rarlee 

n n + i 

linearly with time. The velocity and displacement thus obtained are shown 
in figures 5.15c and 515d. 



indicate the time immediately before and after 

- -f 

the application of the impulse at t . and let v 

f ' 



Suppose t" and t 



n 



and v Indicate 

n n 



respectively the velocity at t~ and t ; these two velocities are related 

by the following equation: 



v n * v n + a n 



The relationship between x , and x . and between x and x^ 

n-l n n n 

are given by: 



X n - X n- l" V n 



Combining equations (5.^7) and (5-^8), the three successive displacements 
are related by: 



This is the basic recurrence formula for the acceleration impulse extrapo- 
lation method. Once the values of x at t , and 



t are known, the 

n * 



, 
n-l 

value at t . can be directly computed without resorting to a trial and 
error procedure. 
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The evaluation of x.. by the recurrence formula needs special con- 
sideration. In equation (5^9)> when n - 0, x.. is given by: 



Both x and x . are equal to zero, and the expression of x. is sim- 
plified to: 



(Atr 



(5-50) 



Hie correct value of a used in the computation is not equal to a 1 ' as 
indicated in figure 5-l6, curve a or curve b, but is given by: 
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Time 
(o) 



At At ' At 1 At 
Time 
(b) 



Figure 5.16. Discontinuities in the acceleration curve 

When equation (5.51) is used in equation (5-50) for the evaluation of X L , 
the value thus obtained is that given by equation (5- 1 * J which is based on 
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the Linear Extrapolation Method. In order to avoid a trial and error pro- 
cedure for the evaluation of x^ , the value of s.^ in equation (5-51) is 

f(t ) f(tj -R^) 

taken to be - instead of its exact value - - - . In figure 

a i 
5.l6a, if a 1 = 0, .a is equal to -5- . In figure 5.l6b if a^ and a 1 

a' 

are approximately equal to each other, the value -=- may be used for a 

Equation (5.^9) is directly applicable for the evaluation of *., x.,... 
x , whenever the acceleration, a(t) , is a continuous curve. When- 
there is a discontinuity in the acceleration this equation is still appli- 
cable provided a modified value of the acceleration is used. This is il- 
lustrated in figures 5-l6a and b. In figure 5-l6a, there is a discontinu- 
ity at ti which occurs at the end of the time interval, At . Under this 
condition, the value of a used in the numerical procedure is the average 
value at discontinuity, namely 2V a h. -+ a k) for \ ' * n fi 6ure 5-l6b, the 
discontinuity occurs within the time interval, At . The correct value for 



a., in this procedure is: 




(5-52) 



For a given problem if there are one or two discontinuities, the 
method described above can be applied conveniently. However, if the number 
of discontinuities is large, it is more convenient to replace the given ac- 
celeration by a smooth curve. 

The main advantage of the acceleration impulse extrapolation method 
is that the numerical computation is straightforward and involves no trial 
and error procedures. Its main disadvantage is that the result is less 
accurate than the linear acceleration method. It has been proved that if 
the acceleration curve is a straight line or a series of straight lines 
with the same slope, the result obtained from the acceleration impulse 
method is exact. However, the actual acceleration curve may be rather 
complicated so that the result obtained by this method is usually only ap- 
proximate. A number of examples have been carried out, and the results 
indicate that the accuracy in the evaluation of the maximum displacement is 
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within 3$ provided the time interval used in the computation is not greater 
than T n /10 . 

Another disadvantage of this method is that any numerical mistake in 
the computation is increased by geometrical progression in subsequent 
steps. Hence, one must be extremely careful to avoid any numerical mistake 
especially in the early stages of computation. 

5-09 NUMERICAL EXAMPIE. An example is given to illustrate the actual ap- 
plication of the different methods of analysis. The same example is eval- 
uated by three different methods: (l) the rigorous method; (2) the linear 
acceleration extrapolation method; and (3 ) the acceleration impulse extrap- 
olation method. 

The external load and the resistance function for the example are 
shown in figures 5-17& and 5-17b respectively. 



f(t) 



o 
o 
o 

.J 




T 
Time 

(a) External Load 



R 



a> 
o 

c 
o 

- 
c/) 



QJ 

a: 



R 



m 



Displacement 
(b) Resistance Function 



Figure 5.17. External load and resistance function for the illustrative example 
The assumed dimensional parameters of the dynamic system are 

P 

m m 2.5 kips-sec /ft 

k M 9860 kips/ft 
R m . 750 kips 
x - 0.0760 ft 

T 0.10 sec 
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B 1000 kips 



2 = 0.10 sec 



The non-dimensional parameters are: 



-0.75, 



Y- = i.oo 

n 



In the rigorous method, the numerical values of the parameters are 
substituted into equation (5.23). The displacement in the time interval, 



< t < t is given by: 



x = 0.01012 



1.012 sin 



H - *'")* H) 



(5-53a) 



The time, t , vhen the displacement x is equal to the maximum elastic 
deflection, x , can be solved from equation (5-53a) and is equal to 0.022 
sec. The displacement and velocity of the mass at t are: 

x(t \ . x = 0.0760 ft 

/ n \ / \ 

01012 (1.012) j~\ cos |2 Y~ - 80 57'J 5-W- ft/sec 
\ n/ \ n / 



= 0. 



Substituting the numerical values of the parameters into equation (5-33*0, 
and integrating, the displacement function in the time interval, 
t < t < t , is given by: 

x = 0.0760 + 5-UU (t - 0.022) + 6(t - 0.022) 2 - 666 (t - 0.022) 3 (5-53*0 
From equation (5-53^), it is found that the maximum displacement is 



x = 0.2825 ft 
m 



which occurs at t given by 

t = 0.077 sec 
m 

The displacement computed from either equation (5-53 a ) or equation (5- 53*0 
is tabulated in column k of table 5-1- 

The same example is evaluated by applying the two methods of numerical 
analysis described in paragraph 5-08. The detailed computations of the 
linear acceleration and acceleration impulse extrapolation methods are 
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given in tables 5-2 and 5-3 respectively. The displacements obtained by 
these three methods are tabulated in table 5.1. It is seen that the error 
in using the numerical methods in the evaluation of the maximum displace- 
ment is less than 1-5^- 



Table 5.1. Comparison of the Results Obtained by 
Different Methods of Analysis 



Time 


Displacement Function x 


t 
(sec) 


Linear 
Acceleration 
Extrapolation 
Method 


Acceleration 
Impulse 
Extrapolation 
Method 


Rigorous 
Method 














0.01 


0.0181 


0.0200 


0.0186 


0.02 


0.0633 


0.0683 


0.0651 


0.03 


0.1182 


0.1216 


0.1203 


0.04 


0.1713 


0.1729 


0.1718 


0.05 


0.2184 


0.2182 


0.2181+ 


0.06 


0.2556 


0.2535 


0.2548 


0.07 


0.2787 


0.27^8 


0.2735 


0-08 


0.2839 


0.2781 










t = 0.077 sec 
m 


0.09 




0.259 1 * 


x = 0.2825 
m 


Maximum 
Displace- 
ment 


0.2839 


0.2781 


0.2825 
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Table 5.2. Details of Computation by the Linear Acceleration Extrapolation Method 





Trial 


t 




V 


- / p t, \/m 
s n ( p n - "n]/- 


Aaeimed 






'= + 1- 


T 


Co^uted 




te. 




(ft) 


(ft/sec) 


(ft/.ec 2 ) 


() 


(kips) 


Ull " J 


U/c ) 


Ul/u.,) 


(ft) 


n 













+1*00 


0.0150 


900 


1*5 


+262 


3.1,1 


0.0180 














o.oiso 


900 


177 


+269 


3*HH 


0.0181 




3rd 










0.0181 


900 


180 


+288 


3-W 


O.OlSl 


i 




0.01 


O.Olfll 


+3. Mi 


+288 


0.0525 


800 


51fl 


+113 


5.w 


0.06*0 




2nd 










0.0640 


BOO 


631 


+68 


5.22 


0.0632 


3 


3rd 


0.02 

0.03 

0.0* 


0.0632 

0.1182 
0.1713 


+5-23 
+5.W 
+5.08 


+70 
-20 
-60 


0.0633 


800 
TOO 
600 
500 


750 
750 
750 


+70 
-20 
-60 
-100 


5.23 

5.M 
+5.06 
+11.26 


0.0633 
0.1182 
0.1713 
0.218>I 


j 




0.05 


0.2l8>t 


+4.28 


-100 




koo 


750 


-iVo 


43.08 


0.2556 







0.06 


0.2556 


+3.03 


-1*0 




300 


750 


-180 


+l.*8 


0.2787 


7 




0.07 


0.2787 


l.*8 


-180 




200 


750 


-220 


-0.52 


0.2639 


the displacement IB larger than the naxlnuB elastic deflection. \ + ^ 1 a constant equal to 750 kips. The 


coBjuitatlnn of f and x j^ beconei itralghtforward vlthout resorting to a trial and error procedure. 


R . kx_ - 9860 x for x < x . 0.0760 ft; R B_ 750 klp f or *w > 0.0760 ft 
n n n c u n 


T n + 1 " 


/"n * "n + 1> 


(AtJ.^.^C* 


/ 


.&)*f 


T n 'V 8 j 


*\3 



Table 5.3. Details of Computation by Acceleration Impulse Extrapolation Method 





t 


P n 


R 

D 


P n-"n 


a . 1? - R \/ 
n V " n " 


. (At) 8 


2x n 


X D _ 1 


X n+l 




(Up.) 


(Xlpo) 


(kijw) 


ft/sec 2 


(ft) 


(ft) 


(ft) 


(ft) 










+1000/2 





+500 


+200 


+O.O2 








tO.O^OO 


1 


+O.01 




+90O 


+197 


+703 


+283 


+0.0283 


+O.OUOO 





+O.06&3 


2 


+O.O2 




+800 


+67!. 


+126 


+50. * 


+0.0050 


+0.1366 


+O.O200 


+0.1216 


3 


+O.O3 




+700 
+600 


+750 
+750 


-50 

-150 


-20 
-60 


-O.O02O 

-0.0060 


+0.2U32 
+0.31*58 


+0.0683 
+0.1216 


+0.1729 
+0.2182 


5 
6 


+0.05 
+O.06 




+500 
+I.OO 


+750 
+750 


-250 
-350 


-100 
-llfO 


-0.01OO 
-0.01>lO 


+0.1.361. 
+O.5070 


+O.1729 
+0.21S2 


+0.2535 
+0.2748 


7 
8 


+O.07 
+0.08 




+300 
+200 


+750 
+750 


-550 


-180 
-22O 


-o.oiflo 

-O.O22O 


+0.51196 
+0.5562 


+6-2535 
+0.27)18 


+0.2781 


9 


+O.O9 






















ta n tor ^n < X e - 


0.0760 


R . 




n 


R - 750 for x > i 

n e 


(At) 0.01 sec 
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5-10 DESIGN CHARTS FOR SIMPLIFIED LOADINGS, a. General. Paragraph 5-10 
presents, in non-dimensional graphs, the results of a systematic analysis of 
single-degree dynamic systems for simplified loads. In order to make these 
graphs simple and convenient for design use, only those cases which involve 
no more than three independent non-dimensional parameters are considered. 
For these cases, the complete results of analysis can be clearly presented 
on a simple graph vith one parameter as the abscissa, a second parameter as 
the ordinate, and the third as a running parameter. 

From the theory of dimensional analysis [l] , the number of independent 
non-dijnensional parameters needed to express the functional relationship of 
all the variables of a given problem is given by equation (5.54). 

N = n - m (5-54) 

where 

N = number of non-dimensional parameters 

n number of dimensional parameters 

m number of fundamental dimensional quantities 

For a single-degree dynamic system, there are three fundamental dimensional 
quantities: namely, mass, length, and time. If consideration is limited 
to problems with three or less non-dimensional parameters, the total number 
of dimensional variables can be no more than six. Of the six variables, 
one is the maximum displacement and one is the mass of the system. Hence, 
the load and the resistance must be completely specified by no more than 
four variables. If two variables, the spring constant, k , and the plastic 
resistance, R , are needed for the resistance function, the dynamic load 
must be completely specified by no more than two variables. These are the 
peak load, B , and the load duration, T '. It is for this reason that 
only very simple types of loads as shown in figure 5-6 are chosen as the 
basis for the preparation of design charts. Although only three simplified 
loads are considered, most actual loads can be replaced by one of these 
with sufficient accuracy for preliminary design purposes. 

The non-dimensional graphs for the linearly elastic, completely 
plastic and elasto-plastic systems subjected to the simplified loads shown 
in figure 5-6 are given separately in the following three sections. The 
non-dimensional parameters chosen for each case are also discussed. For 

39 



5-lCto 
15 Mar 57 

the linearly elastic system, all three types of loads are used but only the 
rectangular and triangular loads are used for the completely plastic and 
elasto-plastic systems. 

t>. Linearly Elastic System. Fora linearly 'elastic system subjected 
to simplified loads shown in figure 5.6, the problem involves five dimen- 
sional variables; x m , m, k, B, and T Hence, only two non- 
dimensional parameters are needed in the presentation of the design charts 
and these parameters are chosen to be: 

x 

D.L.F. Dynamic Load Factor 
x s 



where x o is equal to B/k which is the displAeeaent produced in the 

s 

system when the peak load B is applied statically. 

The method described in paragraph 5-07b can be used for the evalua- 
tion of the dynamic load factor. The displacement function for a given 
load is derived, from which the maximum displacements, x , for different 

values of C_ and the time, t , at which x 's occur, are obtained. 
j. mm 

The dynamic load factors for the rectangular load, the triangular 
load and the step load with rise time, T , are plotted respectively in 
figures 5-18, 5-19, and 5.21. The time ratios t /T are also plotted in 
these figures for each loading. The composite of figures 5-18 and 5-19 is 
given in figure 5.20 for the purpose of comparing the effects of the 
rectangular and triangular loads. 

JC 

The dynamic load factor in this manual is defined as ~ where both 

"* ^J 

x and x are the displacements of the equivalent dynamic system. In 
ni s 

the evaluation of the equivalent dynamic system, it is assumed that the 
deflected shapes of the structure, for example a simply- supported beam, at 
any two instants of time are geometrically similar to each other. Under 
this assumption, the slope and curvature at a given point of the beam de- 
pend only on the mid-span deflection of the beam Thus the dynamic load 
factors also give the ratio of the maximum internal force developed at any 
point under dynamic load to the internal force developed at the same 

kQ 



5-icb m 1110-31*5 -if 15 

15 Mar 57 

point when the peak load, B , is applied statically. 

It should be noted that ordinary structures such as beams or frame 
buildings, involve a continuous distribution of mass. High modes of oscil- 
lation as veil as the fundamental mode which is used in the evaluation of 
the single-degree dynamic system are present. As shown in appendix A, 
EM 1110-3^5-^16, although the higher modes have very little effect as far 
as. the mid-span deflection of beams is concerned, theix effects on the 
bending moment and shear are quite large. Strictly speaking, the dynamic 
load factor based on a single-degree dynamic system, can only be referred 
to as the ratio of displacement but not as the ratio of internal forces. 
For the correct evaluation of internal forces such as bending moment and 
shear, high modes have to be considered and the structure must be replaced 
by a multi-degree dynamic system. In practical application, for the sake 
of simplicity a single-degree equivalent system is still generally used to 
approximate the given structure. The same dynamic load factor is applied 
to the displacement as well as to the bending moment. The high modes of 
oscillation are sometimes partially taken into account in the evaluation 
of the shear and this is discussed in appendix A, EM 1110-3^5-1*16. 

Another point worth mentioning is that the dynamic load factor for 
the equivalent system of a structural member and the dynamic load factor 
for the -quivalent system of a structure as a whole have slightly different 
meanings. For the equivalent system of a member, an elastic beam for 
example, the elastic strain energy is stored all along the beam. If the 
shape of deflection curve remains the same, the internal force at any 
point of the beam is proportional to the mid-span deflection. Hence, the 
dynamic load factor gives the ratio of internal forces under dynamic and 
static loads at any point of the beam. But for the equivalent system of a 
structure as a whole, the strain energy is assumed to be stored in a few 
of its members, for example, the columns of a one-story frame building with 
very stiff girders. For this case, the dynamic load factors give the ratio 
of internal forces under dynamic and static loads at any point of those 
members in which the strain energy of deformation is stored. 

In the application of the dynamic load factor curves of figures 5-l8> 
5.19, and 5.21, the value of C^ is determined from the mass, m , the 
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spring constant, k , and the load duration, T , of the equivalent dynamic 
system. Using C^ , the D.L.F. and t /T are determined from figures 
5.18, 519> and 5-21 for rectangular load, triangular load, and step load 
with rise time, T respectively. 

Tie loading curves determined from the data in EM 1110-3^5-^13 gen- 
erally do not conform to any of the simplified load shapes for which figures 
5.18, 5*19> and 5 21 have been derived. However, with reasonable accuracy, 
these figures may be used for analysis and design where the given load 
curves have been approximated by one of these simplified load curves. The 
method of performing this approximation is explained in paragraph 5-13. 

c. Completely Plastic System. For a completely plastic system sub- 
jected to the simplified loads shown in figure 5-6, the problem involves 
five dimensional variables : x , R , m, B, and T . Only two non- 
dfoensional parameters are needed to express the relationship of- these 
variables in chart form, and these parameters are: 



W 

c w = w 



m 



work done ratio 



C_ = =- = resistance-load ratio 
n JD 



vhere 

W = the maximum work done by the external load and 
m 

W = the fictitious maximum work done defined in equation (5-12) 
The vork done ratios for the rectangular and triangular loads, eval- 
uated by the method described in paragraph 5-0?c, are plotted in figures 
5-22 and 5.23 respectively. The time ratio, t m /T is also plotted in 
these figures for each type of loading. These charts indicate that the 
vork done ratio is independent of the duration of the load, but depends on 
the detail shape of the load; that is, whether the load is rectangular or 
triangular. 

In the application of figures 5-22 and 5-23, the value of W and 
C- are first determined from the peak load, B , the load duration, T , 



the mass, m , and the plastic resistance, R. 



Using C R , C w is de- 



termined from these figures. The maximum displacement, x , is equal 
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* W -/ R where W is the product of C, r . and W 
nr m m y p 

d. Elasto-Plastic System. There are six dimensional variables, x , 
* k > R m B > and T > involved in the problem of an elasto-plastic system 
subjected to a rectangular or a triangular load. Hence, three non- 
dimensional parameters are needed and these are chosen to be 



7or given values- of C T and C R , applying the method described in para- 
graph 5-07d, the displacement function for a given load can be obtained. 
Substituting the maximum displacement into equation (5-38), the maximum 
work done by the external load can be evaluated. This procedure was applied 
to an elasto-plastic system subjected to rectangular and triangular loads. 
The results thus obtained are plotted in figures 5.24 and 5.2?. Because of 
the Importance of the elasto-plastic system in actual application, four 
nore graphs are plotted in figures 5-25, 5.26, 5.28, and 5.29. The first 
two are for rectangular loads and the last two are for triangular loads. 

In figures 5-25 and 5-28, the ordinates are x_/x where x is the 

nr e e 

Buxlmiim elastic deflection. These two graphs are useful in determining the 
Maximum displacement of a given structure for a given load. 

In both figures 5.26 and 5-29, the ordinate is t/T where t is 

nr m 

the time when the maximum displacement is reached. 

In the application of these graphs to an elasto-plastic system, the 
non-dimensional parameters, CL, and C_ are determined from the variables 
of the system. Using C R and C T , the values of CL, , x _/ x and *n/^ 
re determined from figures 5.24, 5.25, and 5.26 respectively for rectan- 
gular load, for figures 5-27, 5.28, and 5.29 for triangular load. The 
Y<mi_iTn work done, W , the maximum displacement, x , and the time, t , 

when x occurs can then be determined, 
m 

Although these graphs are derived for rectangular or triangular 
loads, they are generally applicable to loads which do not conform to 
these simplified shapes. For this condition, it is necessary to idealize 
the given load and this is discussed in paragraph 5-13. 
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figure 5.18. Dynamic load factor and tJT curves for linearly 
elastic system, rectangular load 
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Figure 5.i9. Dynamic load factor and tjt curves for linearly 
elastic system, triangular load 
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Figure 5.20. Composite dynamic load factor and tJT curves for 
linearly elastic system, rectangular and triangular load 
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Figure 5.21. Dynamic load factor and tjl curves for linearly 
elastic system, step load with finite rise time, 7" r 
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Figure 5.22. H'orfe done ratio and t^/T curves for completely 
plastic system, rectangular load 
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Figure 5.23. Work done ratio and tjl curves for 
completely plastic system, triangular load 
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Figure 5.24. Work done ratio for e las to- plastic system, 
rectangular load 
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Figure 5.26> j curves for e las to- plastic system, 
rectangular load 



I 



5 10d 



EM 1110-31*5 

15 Mar 57 



C T = T/T n 



O.I 



o 



Resistonce Function 
Rm 
B 




C T =T/T n 

Figure 5.27. Work done ratio for elasto-plastic system, 
triangular load 
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Figure 5.28. "s m curves for e Las to-plastic system, 
e triangular load 
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Figure 5.23. JLcurves for elasto-plastic system, 
triangular load 
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e. Discussion of the Design Charts. In the preceding paragraphs (b, 
c and d) the response of single-degree dynamic systems subjected to simpli- 
fied loads are presented in non-dimensional graphs. In actual application, 
it would be helpful to knov how the response of dynamic systems varies with 
the peak value, shape and duration of the dynamic load and how the response 
is affected by the mass, spring constant and plastic resistance of the dy- 
namic system. In appendix A, following 'text, a quantitative discussion of 
the effects of small increments of the variables, namely: T, m, k, B, 
and R , on the response are presented in detail and these effects are ex- 
pressed in terms of percentage increment ratios. In this paragraph, only a 
summary of the results derived in appendix A is described. 

First, the response of a linearly elastic system is considered. The' 
dynamic load factor curves for both rectangular and triangular loads are 
shown in figure 5-20. Depending on the value of C^, , these two curvea are 
separated into three regions. In the first region, where C T > 3-0 , the 
D.L.F. for the rectangular load is constant; and the D.L-F. for the tri- 
angular load is practically constant. At a given C^ , the D.L.F. ffinr both 
types of load are approximately equal to each other. Hence, in this region. 
the maximum displacement, x , is: (l) directly proportional to tike peak 
load, B; (2) inversely proportional to the spring constant, k , and 
(3) practically independent of the mass, m , the load duration, T , and 
the detail shape of the initially-peaked load. In the second region where 
C_ < 0.3j the D.L.F. for the rectangular load at a given CL, is aboxrt 
twice the corresponding value for a triangular load. This indicates that 
the area under the load-time curve (the load impulse defined by equation 
(5.5)) is important while the shape of the load is not important. In this 
region, the maximum displacement, x , is: (l) directly proportional to 
the peak load, 3; (2) directly proportional to the load duration, T; 
(3) inversely proportional to the square root of the mass, m; and (k] in- 
versely proportional to the square root of the spring constant, k . In 
the third region, where 0.3 < C m < 3.0, all the variables m, B, T, k, 
and the shape of the load are important. It can therefore be concluded 
that, so far as the maximum displacement of a linearly elastic system sub- 
jected to rectangular or triangular load is concerned, the peak load, B , 
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is alvays an important variable. The spring constant, k , is also impor- 
tant especially when C T is large. The mass, m , and the load duration, 
T , are important only when C T is small. The detail shape of the load 
is important vhen C T is between 0.3 and 3.0. 

Second, the response of a completely plastic system subjected to a 
rectangular or triangular load is considered. The maximum work done, W 
is: (l) inversely proportional to the mass, m; (2) proportional to the 
square of the load duration, T; (3) approximately proportional to the 
square of the peak load, B; and (U) inversely proportional to the plastic 
resistance, R^ . The above conclusions are also true for the maximum dis- 
placement, x , with the exception that the maximum displacement, x 

in 
is inversely proportional to the square of the plastic resistance, R . 

m 
Therefore it can be concluded that all the variables B. T, m and R 

m 
are alvays important in the design of a completely plastic system. 

Third, the response of an e las to -plastic system subjected to a tri- 
angular load as shown in figures 5.27 and 5.28 is considered. Referring 
to figure 5-27* the work done ratio, C^ , is practically constant vhen 
C is smaller than 0.2. In this region, the maximum work done is: pro- 
portional to the square of the peak load] proportional to the square of 
the load duration; and inversely proportional to the mass. Hence, the peak 
load, the load duration and the mass are important variables for short 

duration load. For a given value of CL , the work done ratio is also 

n 

practically a constant for large values of C T when CL is small. Thus 
the peak load, the load duration and the mass are also important variables 
for long duration loads provided CL is small. On the other hand for 

K 

large values of C_ and CL , the values of G W decrease with the in- 
crease of Crp , hence the spring constant is important under these condi- 
tions. The value of C TT also decreases with the increase of CL , hence 

W J 

the plastic resistance, R m , is always important especially at large 
values of CL . 

Similarly, referring to figure 5-28, it can be deduced that, in the 
determination of the maximum displacement, the peak load is by far the most 
important variable. The plastic resistance is important for large values 
of CL provided 0.5 < CL < 1.25- The mass and the load duration are 
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Important for short duration loads and those two variables are also impor- 
tant for long duration loads provided C D is small. The spring constant 

K 

is important for large values of C_ provided C R is also large. 

DESIGN METHODS FOR SINGLE-DEGREE-OF-FREEDOM DYNAMIC SYSTEMS 

5-11 INTRODUCTION. The following paragraphs deal with the design of struc- 
tures or their members wh^ch can be represented by a single-degree dynamic 
system. 

For a given design problem, the dynamic load is always given. Of the 
three variables of the dynamic system, namely, m, k, R , the mass can be 
estimated and hence is considered to be known. There are also the inde- 
pendent variable, time, and the dependent variable, displacement of the dy- 
namic system, the magnitude of which may not exceed specified limits. 
Hence, the designer's problem is to design a single-degree dynamic system 
vith a known mass subjected to a given load such that the maximum deflec- 
tion is equal to or smaller than a certain amount. 

There are two variables, k and R , of the dynamic system which 
muflt be determined to complete the design. Theoretically, there should be 
an infinite number of combinations of k and R to satisfy the given de- 
sign specifications. However, considerations of economy, type of con- 
struction, and standard structural practice limit the relative variation of 
the spring constant, k , and the ultimate resistance, R . Thus, once 
the type of construction has been established, the size or composition of 
the structural members may generally be uniquely determined. 

In paragraph 5-10, non-dimensional design charts are presented for 
idealized loads. The methods of design, using these charts, are first dis- 
cussed in paragraph 5-12. When the given load is different from the 
idealized loadings, the method of approximating the given load so that the 
design charts can be used is discussed in paragraph 5-13. 
5-12 DESIGN METHODS USING CHARTS - IDEALIZED LOADINGS. Three methods 
using the design charts are given in this paragraph for the design of 
single-degree dynamic systems. These are the energy method and the de- 
flection method for cases where plastic deformation is allowed, and the 



5-12* EM 1110-314.5-415 

15 Mar 57 

method based upon the dynamic load factor for cases where only elastic de- 
formation is permitted. 

a. Design for Plastic Deformation - Energy Method. The energy 
method for the design of elasto-plastic systems is based upon the fact that 
the work done on the system up to any time is equal to the sum of the ki- 
netic and strain energies. At the time of maximum displacement, the ki- 
netic energy is zero and the work done by the external load is equal to the 
strain energies. 

When the dynamic load and the mass of the system are given, the fic- 
titious maximum work done on the system by the external load is given by 
equation (5-55)- 

JT 



W 



_ __ 
p ~ 2m 



! 



f(t)dt 



o 
where 



(5-55) 



H - area under the load- time curve 

= BT for rectangular load, or BT/2 for triangular load. 
Suppose, using the method which is described later in this section, 
a trial size of standard structural section is chosen; the value of k and 

R can be determined from the geometrical dimensions of the structure and 
m 

the data in EM 1110-3^5-^1^; then the following non-dimensional parameters 
are computed: 



(5 - 56b) 



For the values of T/T n and C R , the work done ratio, G W = W m /W p is 
then obtained from the design charts. Figure 5-2U is used for rectangular 
load, and figure 5.27 is used for triangular load. The actual maximum work 
done on the system, W ffl , is then computed from the relation W m = (y^ . 
The maximum allowable strain energy to be stored. in the system is 
to the area under the resistance-displacement curve up to the maxi- 
displacement, x ffl . For an elasto-plastic system as shown in 



mum 

59 



BM 1110-3^5-^15 5 " l2a 

15 J'AT 57 



figure 5-7c If the maximum allowable deflection is x ffl , the maximun allow- 
able strain energy storage in the system is 

& ' \ ( x m - X e/ 2 ) (5 ' 57 > 

Comparison of the maximum vork done and the maximum allowable energy 

storage determines the suitability of the^ selected size. 

Based upon the above discussion, a trial and error design procedure 
is set up for an elasto-plastic system as follows: 

Step 1. From the parameters B and T of the given load and the 
estimated mass, the fictitious maximum work done, W , is computed from 
equation (5-55)- 

Step 2. A value of W /W is first assumed and an approximate 

value for W calculated using W from steD 1. A first approximation of 
m p 

the required maximum plastic resistance, R ffl , is determined from equation 
(5.57) using this approximate value for W ffl and neglecting x g as 
follows : 

R m = V X m (5 ' 58a) 

or W +E 

(5-58b) 



where 



m / x \ 

4- - * 



x is maximum allowable deflection, 
m 

Step 3- A trial section is then chosen having a value of R ap- 
proximately equal to the value given by equation (5-58). 

Step U. The actual value of R , the spring constant, k , and the 
maximum elastic deflection, x , are then computed from the selected sec- 
tion and the dimensions of the member. 

Step 5- Factors C and C as expressed in equations (5.56a) and 
(5-56b) are competed. 

Step 6. Using the values of C T and C R , G W is determined from 
the given design chart. Figure 5-2^ is used for rectangular load and 
figure 5-27 is used for triangular load. 

Step 7 The actual maximum external work done, W , is obtained* as 
the product of W (calculated from equation (5.514-)) and C y obtained 
from the chart. 
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Step 8. The maximum allowable strain energy storage in the selected 
trial section is computed by using equation (5-57). 

Step 9- The maximum allowable strain energy storage in the section 
is compared with the actual maximum work done to determine the suitability 
of the selected section. If the allowable strain energy and maximum work 
done are not the same, the above procedure from step 2 to 9 should be re- 
^jeated using equation (5-58b) until a section is obtained which gives the 
closest agreement between these two quantities. In the successive trial 
designs, a new value for G W approximately equal to the value obtained in 
step 6 of the previous trial design is assumed. If the maximum allowable 
strain energy is greater than the maximum work done, the member is over- 
designed, and a smaller structural section should be tried and vice versa. 
The final mass of the system should check the estimated mass in step 1. 

The above procedure has been set up for an elasto -plastic dynamic 
system. The same procedure can also be used for the design of linearly 
elastic systems and completely plastic systems. For a completely plastic 
system, T is equal to zero, hence T/T is always infinite. Thus the 






asymptotic values of C^ for large values of T/T should be used. The 
asymptotic values of C y for figures 5.24 and 5.27 are plotted in figures 
5.22 and 5.23. 

For the design of elastic systems, the curve marked "Elastic" in 
figures 5-2U and 5-27 should be used. In step 2, the approximate value of 
R is computed as 2W /x . The value of CL need not be computed. An 
elastic design is based on dynamic stresses up to the elastic limit of the 
material. If it is desired to limit the stress to less than the elastic 
limit, the computation in step k of R ffl and x g should be based on the 
maximum stress intensity desired in the elastic design. The dynamic "load 
factor method presented in paragraph 5-12c is recommended for elastic de- 
sign. This method is based on the computation of equivalent static loads 
or deflections as opposed to energy considerations. 

The above procedure gives an exact solution where the external load 
is either rectangular or triangular, because the charts from which C^ 
are obtained are computed based on these two types of loads. For any 
other type of loading which can be completely specified by two parameters, 
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B and T , similar charts can be prepared. The above procedure would 
still give the exact solution provided the value of C^ were obtained from 
the design chart corresponding to the given type of loading. However, in 
actual problems, the loading is rather complicated and generally cannot be 
specified by two parameters. It is impractical to prepare design charts 
which will give the exact solution for any type of loading. Hence, the de- 
sign charts are only prepared for two basic types of loading. For an actual 
problem, the loading must be approximated by either a rectangular or a tri- 
angular shape. The method of approximating the given load so that the de- 
sign chart for rectangular or triangular load can be used is discussed in 
paragraph 5-13- 

The application of the energy method of design to actual problems is 
given In EM 1110-3l*-5-4l6 and kYf. In paragraph 6-11, EM lllO-S^-lnS, a 
numerical example is given to illustrate the design of beams. 

b. Design for Plastic Deformation - Deflection Method. The deflec- 
tion method of design is practically the same as the energy method with the 
exception that the non-dimensional parameter, x /x , is used instead of 

the vork done ratio, C . The design procedure for the deflection method 

w 

is as follows: 



Step 1. A value of C_ ^..B is assumed, and a first approximate 



value of the resistance, R , determined. 

7 m 

R m = C R B (5.59) 

Step 2. A trial section is chosen having a plastic resistance, R , 
approximately equal to the value given by equation (5-59). 

Step 3. The actual values of R m , the spring constant, >k , and the 
maximum elastic deflection, x , are computed from the selected section 
and the dimensions of the member. The mass, m , of the system is also 
determined. 

Step U. Factors C T and C R as expressed in equations (5-56a) and 
(5-56b) are computed. 

Step 5- x ji/ x e is determined from the given design charts. Figure 
5.25 is used for rectangular loads, and figure 5.28 is used for triangular 
loads. 
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Step 6. The maximum deflection, x m , is obtained as the product of 
** x / x ' 



Step 7- Comparison of the computed maximum deflection, x , with 

the design specification for the maximum allowable x determines the 

m 
suitability of the selected trial section. By assuming a new value of C 

the design procedure from steps 1 to 7 is repeated until a section which 

gives the closest agreement between the computed x and the allowable 

m 

x is obtained. 

is 

c . Design for Elastic Deformation - Dynamic Load Factor Method . 
This design procedure is for the case in which the external load can be ap- 
proximated by either a rectangular or triangular shape or a step load with 
a finite rise time as covered by the design charts in figures 5.18, 5.19, 
and 5.21. Methods for approximating a given load to conform to one of 
these simple types are covered in paragraph 5-13- 

Step 1. A value of dynamic load factor is first assumed. 

Step 2. A trial section is selected having a maximum resistance, 
R , approximately equal to 

R m (D.L.F.)B 
where 

D.L.P. is the assumed dynamic load factor. 

Step 3. The spring constant, k , is determined as described in 
paragraph b above, and the maximum elastic deflection is determined for 
various type elements as described in paragraphs 6-13 through 6-2^ in 
EM 1110-3^5-^16. _ 

Step U. The time ratio, C T T/T n = g Vm is determined - 

Step 5. Using C , a revised value for the dynamic load factor is 
determined from figures 5.18 or 5-19 for rectangular or triangular loads 
respectively, or from figure 5.21 for step loads with a finite rise time. 

Step 6. The static deflection, x g , is determined by x g = B/k . 

Step 7. The maximum deflection, x m , is obtained as the product of 
x and the revised dynamic load factor. 

Step 8. Comparison of x ffl and x g determines the suitability of 
the selected trial section. The design procedure from steps 1 to 8 is re- 
peated until a section which gives the closest agreement between x m 
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5-13 



and x is obtained. If it is desired to limit the maximum deflection to 
less than the maximum elastic deflection, a fictitious x g equal to the 
desired maximum deflection can be used. The maximum resistance of the sec- 
tion, R , would be reduced in the same ratio as x to correspond with 
in 

the smaller design deflection. 

An alternate dynamic load factor method may be used if desired, in 
which stresses are used rather than deflections. To use this approach, 
proceed as above but omit computation of x , x , and x . For step 8 

compare R of section with (D.L.F.)B to determine suitability of the 

m 

selected trial section. By either approach, it is necessary to determine 

t and check the load shape as discussed in paragraph 5-13- 
m 

The application of this method to the design of elastic beams and 
one-story buildings is illustrated by numerical examples in EM 1110-3^5-14-17. 
5-13 METHOD OF APPROXIMATING A GIVEN LOAD. The design methods given in 
the previous paragraph are based upon an idealized load acting on a dynamic 
system with an idealized resistance function. For an actual problem, the 
resistance function is very close to one of the idealized forms for an 
elastic or an elasto-plastic system; however, the given load may be radi- 
cally different from either a rectangular load or a triangular load or a 
step load with a finite rise time. Hence, the first step in the design is 
the replacement of the given load by one of the idealized loads for which 
the design charts are provided. 

The basic consideration in replacing the given loading by an ideal- 
ized loading shape is that the maximum displacement produced by both must 
be equal. This replacement of load usually introduces a certain amount of 
error and is one of the more difficult parts of the whole design problem. 
The amount of error depends largely upon how well the loading is approxi- 
mated which depends entirely upon the judgement of the designer. No gen- 
eral rules can be given for replacing a given load by a simplified load; 
however, the following is set up as a guide. 

Suppose for a given problem, the load is shown in figure 5-30a. 
In the first approximation a rectangular load as shown by the dotted curve 
may be used. Following the design procedures set up in paragraph 5-12, the 
values of C T and C R are obtained. Then the value of t /T is obtained 
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from figure 5.18, 5-22, or 5.26. This parameter is very important in 
Judging how good the load approximation is. Since t is the time at 



m 



which x m occurs, the approximate load and the given load should (l) be 
very close to each other in the time interval of time from to t , and 
(2) have about the same area in this time interval. 

Depending upon the value of t obtained from the first approxima- 
tion, the given load should be replaced by one of the rectangular loads 
shown in figure 5-30b and figure 5-30c. However, if t ffi is much larger 
than Tp as shown in figure 5.30a> it is expected that the internal re- 
sisting force in the time interval from to T is small. Based on the 
discussion in paragraph 5-05a, for this case, the detail shape of the load 
is not important, and the given load can be replaced by a pure pulse load 
with equal total impulse. In this case, the value of C is approximately 

W 



to unity, hence W ffl = W 
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Figure 5.30. Typical load approximations 
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In figure 5-30a, b, c, a rectangular load is used to replace the 

given load, hence the value of t is obtained from figures 5-18, 5-22, 

m 

or 5.26. For the case when the given load is replaced by a triangular load 
(figures 5.30d and 5.30e), t may be obtained from figures 5.19, 5.23, 
or 5.29. When the given load is replaced by a load with a finite rise 
tine (figure 5.30f ) t may be obtained from figure 5-21. 
5-llf DISCUSSION OF DESIGN METHODS, a. Accuracy of Design Methods. From 
the above discussion, it is seen that in applying the design procedures 
given in paragraph 5-12 to an actual problem, idealized load and resistance 
values have to be used to replace the actual values. It is expected that 
a certain amount of error will be involved. A final check of the actual 
design using the numerical methods of analysis described in paragraph 5-08 
is necessary. Therefore, this design method can be considered as a pre- 
liminary design procedure. The result obtained from the preliminary design 
must be checked by a numerical method of analysis using the actual load and 
resistance function. 

The maximum displacement evaluated by the numerical method may be 
different from the design specification. There are two reasons for this 
discrepancy. One is the fact that idealized loads and resistances are used 
in the preliminary design. The error from this source is usually less than 
5/&. The second is due to the use of standard structural sections. Only a 
limited number of sections are available. The difference in properties be- 
tween two adjacent sections may be more than 10$. For a given design prob- 
lem, it is sometimes found that no standard section is available in order 
to satisfy the given design specifications. If a standard section close 
to the desired section is used, the maximum displacement may be different 
from the design specifications by more than 10$. However, a number of de- 
sign problems have been carried out and it has been found that the accu- 
racy of the preliminary design method, judging by the difference between 
the actual and the specified maximum displacements, is within 10$ for the 
design of one-story buildings, and 25$ for the design of structural mem- 
bers, such as beams or slabs. The larger error in the design of beams and 
slabs is due to their shorter periods which make them more sensitive to the 
shape of the load, and to gaps in available structural sections (see 
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appendix A following the text of this manual). 

When designing systems with an elasto-plastic -behavior, a small change 
in the maximum resistance of the section will cause proportionately larger 
changes in the maximum displacement (see appendix A). Because of this 
sensitivity, a 100 error in the displacement represents a relatively 
smaller change in the required maximum resistance. Therefore, numerical 
analysis is not always justified in checking the design of an elasto- 
plastic system. 

b. Comparison of Different Methods of Design. In paragraph 5-12, 
three different methods of design are given. The energy and deflection 
methods. are applicable to the design for plastic deformation while the dy- 
namic load factor method is applicable only to the design of elastic 
systems . 

In both the deflection and the dynamic load factor methods, the non- 

dimensional parameter used is a ratio of deflections, x /x , for the 

m' e ' 
former, and x m /x s for the latter. Hence, the dynamic load factor method 

is actually only a special case of the deflection method. 

Basically, there are two different methods of design, the energy 
method, and the deflection method. The design procedures for these two 
methods as given in paragraph 5-12 are on a trial and error basis. After 
a trial section is selected, the amount of numerical work involved in de- 
termining the suitability of the selected section is about the same no 
matter which method is used. But the convergence of successive trials for 
a given problem depends on which method is used. 

In the design of one-story buildings which can be replaced by a 
single-degree dynamic system, the value of C_ is of the order of unity; 

and the value of C_, varies from 0.25 to 1.5 for an elasto-plastic system 
JK 

under triangular load. In these ranges of values for CL and C^, , the 
maximum work done, W , is not very sensitive to small variations in the 
spring constant, k , or the plastic resistance, R m (appendix A). In 
other words, the maximum work done is practically constant when different 
sections are used in a trial and error design procedure. Thus, in suc- 
cessive trials, the only modification is in the strain energy storage. If 
in one trial design the section is found to be unsatisfactory, a second 
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trial can be started by using another section- with maximum strain energy 
storage approximately e^ual to the work done obtained in the previous trial 
design. On the other hand, if the deflection method is used within the 
same ranges of values of C R and T/T R as mentioned before, the maximum 
displacement is rather sensitive to any change in the value of the spring 
constant, k , or plastic resistance, R m . If a trial design is found to 
be unsatisfactory, there is practically no guide as to what the next trial 
section should be. Hence, in the design of one-story buildings, the energy 
method of preliminary design is preferred to the deflection method. 

However, in the design of individual members such as beams or slabs 
under triangular load, the value of C T is usually greater than 3 and, the 
value of C is usually greater than 1. In this range of values of C T 
and C , the variation of W ffl with any change in k or R ffl is about of 
the same order of magnitude as the variation of x ffl (appendix A of this 
manual). Hence, in the design of beams or slabs, the deflection method is 
as convenient as the energy method so far as the convergence of successive 
trials is concerned. 

So far the discussion has been limited to the design of a single- 
degree dynamic system. For a mult i -degree dynamic system such as a multi- 
story building, the design problem is complicated by the inter-actions be- 
tween different stories. It is shown in paragraph 5-21c that, for design 
purposes, any story of a multi-story building can be replaced by an equiv- 
alent single-degree dynamic system. From this, a method of design based on 
energy considerations has been developed which gives reasonably good 
results. 

ANALYSIS OF MULTI-DEGREE-OF-FREEDOM DYNAMIC SYSTEMS 

5-15 INTRODUCTION. In the following paragraphs, the general problem of 
analysis of mult i -degree -of -freedom dynamic systems is considered. The 
types of structures and structural members which are replaced by multi- 
degree dynamic systems for analysis are, for example: (l) beams under 
lateral load when higher modes as well as the fundamental mode are consid- 
ered (appendix A of EM 1110-31+5-1*16); (2) the effect of support movement on 
the response of a beam; and (3) multi-degree buildings under lateral loads. 
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In order to make the problems under discussion easier to be visual- 
ized, the terminology of multi-story buildings is used in this and the sub- 
sequent sections of this manual. It should be understood that the methods 
of analysis are applicable for any multi-degree dynamic system and are by 
no means limited to multi-story buildings. In paragraph 5-16, the equa- 
tions of motion are set up. The rigorous and numerical methods of analysis 
are given in paragraphs ^-17 and 5-18 respectively. 

5-16 EQUATIONS OF .MOTION. For the dynamic analysis of a multi-story build- 
ing of n stories acted upon by lateral blast loads, it is necessary, to re- 
place the structure by an equivalent dynamic system. The distributed masses 
of the structure are lumped together at the floor and roof levels, giving a 
series of n lumped masses. The masses are connected to each other by 
weightless springs which simulate the lateral resistance of the structure. 
Similarly the distributed loads on the structure are lumped together at the 
lumped masses. Hence, the multi-story building of n stories is repre- 
sented by an n-th degree dynamic system with n concentrated masses con- 
nected by weightless springs and acted upon by concentrated loads at each 
of the masses. The variables under consideration are the displacements of 
these 'masses. 

The following notations are used: . . 

For masses: m^, nig....m _ ^, m , m + ^....ra 

For absolute displacement of these masses: 

X l' V"' X g - 1' V X g + l"" X n 
For external load on these masses: 

f^t), f 2 (t)....f g _ x (t), f g (t), f & + I (t)....f n (t) 
The equation of motion of mass, m , is given by equation (5-60): 

O 

d 2 x 
"g-^-V** - R gt (5 ' 60) 

g at s K 

where R ,. is the total resisting force which is exerted by the columns 

gt 
and walls above and below upon the mass m . Depending upon the type of 

6 
construction, two cases are considered for the determination of R^ . For 

example, in the case of a shear wall structure, where the building is con- 
sidered as a vertical cantilever and the effects of over-all bending on the 
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building is taken into consideration, the resistance function, R^ , de- 
pends upon the displacements of all the n masses. In the case of frame 
structures where the effect of over-all bending of the building is neglected, 

R is a function of only three absolute displacements, x , x , and 

gt & - J- g 

a. Multi-Story Shear Vail Buildings. The resisting forces acting on 
any floor of a multi-story shear wall building depend on the displacements 
of all the floors. Figure 5-31 shows the absolute displacements of any two 

floors represented by two masses, 
m and m. of a multi-story build- 
ing. The absolute displacement of 
^ _ each mass is represented by x and 
x. respectively, their relative dis- 




placement being given as 

Due to the difference in magnitudes 



x - x 

& X 



of the displacements x. and x , 
a pair of internal forces, R . , as 
shown in figure 5-31 are developed. 

The masses m. and m are said to 

i g 

be coupled together. If the coupling 



spring constant is k 



. 



the expres- 



sion for the resistance function, 



R 



in the elastic range is given 



by equation (5-6l). 



Figure 5.31. Any ttuo floors of a 
multi-story building 
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The positive direction for R 
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gi " gi 
as well as the positive directions for 



(5-61) 



x are indicated by arrows in figure 5-31. 

O 

For & shear wall building, the coupling of each floor mass to every 
other floor mass and the ground must be considered. The general expression 
for the total resisting force, R , acting on the g-th floor mass of a 

6& 

mul^i-story shear wall building is: 



TO 
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The design of multi-story shear wall buildings will generally be simplified 
to elastic analysis because of the small allowable lateral deflections anG 
the complexity of the problem. For this case, the above expression may be 
written as: 



- + - x * " * - * + X - *) 






In equation (5.62), R is the resistance developed by the shear wall at 
the g-th floor due to the relative movement of m with respect to the 
ground. Note that when i = g, R . = R = . In equation (5.63), the 
ground displacement is assumed to be zero. 

Substituting equation (5-62) into equation (5.60), the equation of 



nction for mass, m , is: 

O 



.2 
d x 



m _ 1 = f (O - > K (5-6U) 

s dt g /-' gl 

1=0 
In the elastic range, equation (5. 6*0 becomes 



i = 

There is one equation of the same type as equation(5-6 1 *) for each floor 
mass (see figure 5-32). Hence, the equations of motion for an n-story 
building consist of n simultaneous second order differential equations. 
* The resistance function, R . , may also be expressed by: 

n 

T? v'v j-ir'v j. v'v K^ V' v * wh^re k' . ~*^ * and. 
gt ~ gl 1 g2 2 + K gn x n " L, ^ K gi i' gi 8 1 

k gg * k go + k gl + V + k g3* ' ' ' k gi + " " k gn ' 
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Figure 5.32. Equivalent dynamic system for an n-story shear wall building 

For a multi-story shear vail building in the elastic range, the 
variables of the equivalent dynamic system can be evaluated based on methods 
which are given in detail in EM 1110-3^5-^16 and EM 1110-3^5-^19. 

Having determined the variables of the system, and the load acting on 
the structure, the analysis is carried out for the elastic range using 
equations of the form 5-65, as explained in paragraphs 5-1? and 5-lB. 

As soon as the stress at any point of a shear wall for any story ex- 
ceeds the elastic limit, not only the internal resisting force in this 
story is affected, but also practically all of the coupling constants of the 
equivalent dynamic system are affected. Whenever the stress in any story 
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exceeds the elastic limit, the resistance functions of the complete struc- 
ture have to be re-computed. Hence, it would not be practical to attempt 
to analyze a multi-story shear wall building with stresses exceeding the 
elastic limit. The application of the above equation is therefore Halted 
to the analysis of shear wall buildings in the elastic range. 

b. Multi-Story Frame Buildings. 
In frame buildings, most of the re- | 

sistances, R . , in equation (5-62) 
are either zero or negligibly small 
except for the resistances, R 
and R 

ances developed between adjacent 
iMLSses. In an actual building, 
these resistances are developed be- 
tween the floors by the columns and 
walls. Because of the relatively 
large stiffness of the frame ele- 
ments in the plane of the floor in 
frame buildings, the lateral re- 
straints are not so important be- 
tween non-adjacent floor levels. 
Eliminating all resistances except 
those developed between adjacent 
masses, the equivalent system in 
figure 5.33 is obtained. This 
equivalent system, which is appli- 
cable for both elastic and elasto-plastic behavior, is used in the anal- 
ysis of multi-story frame buildings. 

For the equivalent system shown in figure 5-33> using the following 
simplified notations: 




n,n-l = Up 
i 

n,n-l 

R n-l,n-2 = R n -| 
R g+2,g-H = Rg+2 

R g -t-l.g - Rg +| 



R I0 5R I 



Figure 5.33. Equivalent dynamic system for 
o multi-story frame building 



R 



g - r 



the total resisting force, R , , on 



"lO 



mass, 



g, 



m 



is 









R g - 



(5-66) 
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In the elastic range, equation (5-66) becomes 



O \r Iv v 
R gt = g( X g " g - 

Substituting equation (5-66) into equation (5-60), the equation of 
motion for mass, m , is given by 

o 

d x e 
m g f = f g (t) - R g + R g + I (5-68) 

In the elastic range, equation (5-68) becomes 
^2 



Rearranging terms in equation (5-69)* tne following recurrence formula 
expressing the relationship of any three consecutive displacement variables 
is obtained. This is given by 



d 2 x 



( k g + k s + i) x g k s + i x g + 1 - f g (t) (5 - 70) 



For an n-th degree dynamic system as shown in figure 5-33* there are 
n separate equations of the same type as equation (5'70). These n equa- 
tions are 

AI 



at 



/, , \ 

Ik 4-Kiv -Vv _ V v -P i -t- i 

i^2 T ^j 2 21 ^^o " I 2\ T '^ 

: ' * ' (5 . 71) 
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d 2 x 






n 






ffi. '" JT " 

n dt 2 


n x n " k n x n - 1 = 


f n (t) 



In equation (5-71), all the coupling springs are assumed to be in the 
elastic range. If a certain spring such as k g has reached its plastic 

resistance, R^ _, this value should be used for R instead of 

6 



. 

The analysis of an n-th degree dynamic system involves the evaluation 
the n simultaneous differential equations. 

5-17 RIGOROUS METHOD OF ANALYSIS. In the rigorous method, the standard 
procedures of solving simultaneous differential equations are used. Load 
and resistance functions must be expressed in simple mathematical forms. 
The displacement of any mass, say x , is first evaluated. The remaining 
x's are then determined using the recurrence formula given by equation 
(5.70) for frame buildings, or equation (5.65) for shear wall buildings. 

Eliminating all the x's except x from the simultaneous equations 

D 

given by equation (5-71) for frame buildings or equation (5.65) for shear 

vail buildings, a linear differential equation of the (2n)th order in x 

g 
is obtained. Just like the case of single-degree dynamic systems, the 

solution of x of an n-th degree dynamic system consists of two parts. 

O 

One is the transient term and the other is the forced solution. Thus the 



complete solution of x can be expressed as 

O 



x 



g 



(5-72) 



where x is the total displacement of the mass, m ; x ,, is the portion 

O ' O O 

of x which is referred to as the forced solution; and the term in 

O 

brackets is the portion of x referred to as the transient solution. 

6 

There are 2n constants of integration in the complete solution. These 
constants are determined from the initial conditions that all the n massei 
are at rest at the initiation of the load. 

There are n oscillatory terms in the transient solution. The fre- 
quencies of oscillation are the n natural frequencies of the n-th degree 
dynamic system. The lowest frequency is defined as the fundamental fre- 
quency of the dynamic system. 
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Several other methods, based on the same principle as described 
above, but differing in detailed techniques are also commonly used in the 
analysis of multi-degree dynamic systems. In one method, the set of simul- 
taneous differential equations (5-71) or (5-65) is first transformed into a 
set of algebraic equations by applying the principle of laplace Transforma- 
tion and operational calculus [^]. The complete solution is then obtained 
by inverse transformation. In another method, the solution of the problem 
is started from the recurrence formula, equation (5-70), by applying the 
methods of solving finite difference equations [8j. 

There is another method which is radically different from the above 
methods. In this case, the principle of generalized coordinates is used 
[2, 6, 7]. An n-th degree dynamic system is changed to an equivalent system 
with n uncoupled masses. Each mass is individually supported by a spring 
and acted upon by a generalized load. Thus the n-th degree dynamic system 
is changed to a system of n uncoupled single-degree dynamic systems. 
This method of generalized coordinates is extensively used in the study of 
structural members such as beams, columns, or slabs, when the high modes 
of oscillation are to be considered. The application of this method to the 
analysis of elastic beams is given in appendix A, EM 1110-3^5 -Ul6, and to 
the analysis of multi-story frame buildings in reference 6. 

The difficulty in applying the rigorous method of analysis to multi- 
degree dynamic systems is evident. The problem involves the solution of n 
simultaneous differential equations of second order with 2n initial con- 
ditions. Suppose two separate equations are needed for each load and the 
dynamic system is elasto -plastic; then, as many as ^n sets of equations 
are needed to express the equations of motion. In each set, there are n 
simultaneous equations. Thus, in general, the rigorous method of analysis 
is impractical for multi-degree dynamic systems, and numerical methods must 
be used. The two numerical methods described in paragraph 5-08 for single- 
degree dynamic systems are also applicable to multi-degree systems. These 
two methods are described separately in paragraph 5-18. 

5-18 NUMERICAL METHODS OF ANALYSIS, a. Linear Acceleration Extrapolation 
Method. The equation of motion for mass, m , is given by equation (5-6o). 

O 

Suppose a denotes the acceleration of m ; the expression for a as a 
function of time is 
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d 2 x g , 

&~ = ~ O = "~~~ If It) -R (til f C T5 \ 

g dt 2 m g Le . gt v J ] t5-73J 

where ^(t) is the total internal resisting force acting upon m at 'any 
time, t . The expressions of R for the shear wall building and frame 
building discussed in paragraph 5-16 are given respectively by the follow- 
ing equations: 



n 



Shear wall building, R t = R t (5.7*0 



Frame building, R gt = R g - R g + ^ ' (5.75) 

If the acceleration a is assumed to be linear in the time interval 

D 

from t n to t n + l , from equation (5-410, the relationship between the 
displacements of x at t R and t + , is given by 

- + (At) 



where /dx\ 

t n = \ aT/t = t 

At - * ' * 



( v s) 



Based upon equation (5-?6), a trial and error procedure can be set up for 
the computation of x's at t from the values of x's at t . A set 

of n values of x, , x_, ---- x . ---- x are assumed for the displacement of 
1* 2 g n 

each mass at time t , from which the resistance R's and the accelera- 

n + 1 

tion a's at time t , are determined. The actual values of x's at 
n + 1 

*n + 1 ^^ then com P uted from equation (5-76). Comparison of the computed 
and assumed values of x's indicates the correctness of the assumed values. 
This trial and error procedure is tedious to carry out since it 
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involves n variables. The computed values of all the n variables must 
agree respectively with the n assumed values. On the other hand, this 
method usually gives accurate results. When all ( the coupling springs are 
in the plastic range where the resistances are constant, the computation 
becomes straightforward. Further details of this method are given in para- 
graph 5-08 where it is described for the analysis of a single-degree dy- 
namic system. 

b. Acceleration Impulse Extrapolation Method. The recurrence for- 
mula, equation (5. ^9), for the single-degree dynamic system is also appli- 
cable for multi-degree dynamic systems. For the displacement, x , of 

o 

'mass, m . the recurrence formula is 
g 



t a 



where 



At = t , - t = t - t , = = t_ - t, 

n + 1 n n n - 1 2 1 

(a V _ = the acceleration of mass, nig , at time, t = t , 
/ ~ n the expression given by equation (573). n 

The method of modifying the acceleration at any point of discontinuity 
of the acceleration curve is the same as in the case of single-degree dy- 
namic systems (see paragraph 5-08). The advantage and disadvantage of this 
method as described for the single-degree dynamic system also hold for 
multi-degree systems. If the time interval (At) is chosen to be about 
T Q /10 where T R is the shortest period, of oscillation, the accuracy in 
the maximum absolute displacement is about 3%. 

The application of the acceleration impulse extrapolation method to 
the analysis of multi-story frame buildings is shown in EM 1110-3^5-^18. 
The equivalent dynamic system for frame buildings is the case shown in 
figure 5-33- 

DESIGN OF MULTI- STORY FRAME BUILDINGS FOR PLASTIC DEFORMATION 

5-19 INTRODUCTION. The problem considered in this and the following para- 
graph is the design of a multi-story frame building which is represented by 

78 



'5-19 Bi 1110-345-415 

15 Mar 57 

the equivalent dynamic system shown in figure 5-33- The terminologies of 
multi-story buildings are used in the discussion. Thus the concentrated 
masses are approximately the masses of different floors, the displacements 
O f each mass are -the displacement of each floor level, and the springs are 
the columns between adjacent floors. 

In the design problem, the floor masses can be estimated and the ex- 
ternal loads acting on the masses are given. The maximum allowable dis- 
placements between adjacent floors are prescribed. The problem is the se- 
lection of the column sections for each story such that the relative dis- 
placements between adjacent floors do not exceed the prescribed maximum al- 
lowable values. 

One approach to the design problem is the preparation of non- 
dimensional design charts as in the case of single-degree dynamic systems. 
Once charts are available, the complete design can be carried out in a 
series of routine computations. A little reflection will indicate that, 
even for a two-story building, there are too many non-dimensional param- 
eters to permit the preparation of practical design charts to cover all pos- 
sible cases. Thus the design of multi-story buildings cannot be carried 
out in as simple and straightforward a manner as the design of single- 
degree dynamic systems. 

Another method of approach to the design problem is to first assume 
the sizes of all the columns for the entire building and then check the 
adequacy of the assumed sizes by a numerical analysis of the whole struc- 
ture. Any readjustment of the column size is then made according to the 
results of the numerical analysis. The size adjustment is difficult be- 
cause the change of column size of any one-story affects the relative dis- 
placements of the entire structure. If, in the numerical analysis, the 
relative displacements between several adjacent floors are found to be un- 
satisfactory, it is difficult to decide which column needs modification, 
and by what amount the size of the column should be changed. The detailed 
numerical analysis may have to be carried out over and over again before a 
set of satisfactory column sizes is found. This method, therefore, is 
difficult and tedious to carry out, and is not recommended except for 
problems where a rational method of design has not been developed. 
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The design method recommended for multi-story frame buildings is baaed 
on a floor-by-floor iteration procedure in which the column size for the 
first, second, third, and on up to the top story is selected one after the 
other. Specifically, the outline of the design procedure is as follows: 

Step 1. The plastic resistance and the size of the columns for the 
first story are determined on the assumption that (l) the first floor is 
a single -degree dynamic system with a mass, m^ , subjected to an external 
load, f-^t) , and (2) all the remaining floors are absent. 

Step 2. The plastic resistance, R^ , of the columns for the second 
story is determined by (l) incorporating the result of the preliminary de- 
sign of the first story, and (2) assuming that all those floors above the 
second floor are absent. 

Step 3. The plastic resistance, R,. , of the columns for- the third 
story is determined by (l) incorporating the results of the preliminary de- 
sign of the first and second stories, and (2) assuming that all those 
floors above the third floor are absent. 

Step 4. The above procedure is continued until the ultimate resist- 
ance, R , of the columns for the top story is determined. 

Step 5. The preliminary plastic resistances determined in step 1 
through k for each story are adjusted by taking into account the effect on 
the columns of a given story by the columns of the story which are directly 
above the given story. For example, the shearing force developed in the 
columns of the third story is considered in modifying the plastic resist- 
ance of the second story. 

Step 6. The column sizes for different stories are then determined. 
The sizes of the selected columns should be based on the revised resistance 
values obtained in step 5. In the determination of the size of the se- 
lected columns, the effect of vertical load on the plastic moment of the 
section should be taken into account. 

Step 7- A numerical analysis for the entire building is then car- 
ried out to determine the suitability of the selected columns. Any revi- 
sion of the column size can be made according to the result of the numer- 
ical analysis. 

There are two basic procedures in the above floor -by-floor iteration 
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design method. One is the preliminary determination of the resistance of 
the columns for any story after the column resistances of all those stories 
below this story have been determined. The other is the adjustment of the 
preliminary resistance thus determined for any story after the column re- 
sistances of all those stories above have been determined. The fundamental 
principles in carrying out these two basic procedures are discussed in 
paragraphs 5-21 and 5-22 respectively. The complete design procedure for 
the design of multi-story buildings is given in detail in paragraph 5-23. 
Die assumptions used in the floor -by-floor design procedure and the ac- 
curacy of this design method are summarized in paragraph 5-25. 

In paragraph 5-20, the absolute displacements, relative displace- 
ments and resistance functions of a three-story building are plotted versus 
time- These curves are used (l) to illustrate the general characteristics 
of multi-story buildings, and (2) to explain why a somewhat complicated 
procedure has to be used in design. 

5-20 DISPLACEMENTS AM) RESISTANCE FUNCTIONS OF A THREE- STORY FRAME 
BUILDING. A three-story frame building designed for blast resistance is 
illustrated in figure 5- 3^- The floor-by-floor resistances and displace- 
ments of this building under the design blast loading are plotted in 
figure 5.35. This building was designed based on the specification that 
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Figure 5.35. External loads, displacements and resistances of three-story 
frame building as functions of Lime 

the maximum allowable relative displacement for any story should not exceed 
0.33 ft- Subscripts 1, 2, and 3 refer to the first, second, and third 
floors above the ground respectively. Hence, the third "floor" is the 
roof. 

Due to the differences in the floor masses and in the external loads 
on the masses, the relative displacements between m, and ni , >' = 
X 2 " X l ' and tetween m o and - m o > x o - x o " x o > are negative shortly 
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after the initiation of the load. The relative displacement, X , is zero 
at t and the relative displacement, ^ , 8 zero at t 1 . In the time 
interval from to t , the resistance, ^ , is negative; and in the time 
interval from to t^ , E^ is negative. The maximum elastic deflections 

for the first, second, and third stories are reached at time t t 

le ' 2e ' 
and t respectively. 

The curves in figure 5-35 show three important characteristics vhich 
are generally true for any multi-story frame building. First, the ratio of 
the absolute displacements of any two floor masses varies over a rather 
large range. Hence, in design, it cannot be assumed that the absolute dis- 
placements of any two floor masses always bear a constant ratio. Secondly 
the resistances, Rg , or R^ are negative in a short interval of time 
after the initiation of load. The subsequent displacement of nu or nu 
is considerably affected by the negative portion of the resistances. Hence, 
in design, it is necessary to take the negative portion of the resistance 

function into account. Thirdly, the maximum absolute displacements, x 

1m 
x , and x.. of m.^ nig, and m_ do not occur at the same time. As 

shown in figure 5-35b> x^, x^, and x^ occur at t^, t^, and t 
respectively, and t < t < t . When m^ reaches its maximum dis- 
placement, x , the floor mass, m, , is moving backward due to elastic 
rebound. Hence, in the evaluation of the maximum relative displacement be- 
tween the first and second floors, it is necessary to consider the elastic 
rebound of floor mass, EL . 

The floor-by-floor iteration design procedure presented in these para- 
graphs is Intended to take these important characteristics into account. 
In the preliminary design of the columns for any story, the resistance as 
a function of time for the story must be estimated. The energy method of 
design given in paragraph 5-12 for single-degree dynamic systems, with 
slight modification can be applied to the determination of the required 
plastic resistance for the story. In the evaluation of the maximum allow- 
able strain energies of the columns, the maximum absolute displacements for 
masses, m-,, m 2 , and m_ given in equation (5.?8) niust be used. 

(5 - 78l) 
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vhere X 



, , 
1m 



, and X_ are the prescribed maximum relative displace- 
2m 3 m 



merits for the first, second, and third stories respectively] and X.^ , X 
are respectively the maximum elastic relative deflections for the first and 
second story columns. 

In equation (5.?8), the maximum absolute displacements x.^, x^, 

x, are expressed in terms of the maximum allowable relative displacements 
3m 

by assuming that the rebound of the lover stories equals their maximum 
elastic relative deflection. The maximum absolute displacements thus de- 
termined are used for the computation of the energy absorption by equations 
given in paragraph 5-21c. 

5-21 THE PRELIMINARY DESIGN OF CQUJMNS FOR ANY STOHY. a. Effective Load 
and the Equations of Motion in Terms of Relative Displacements. In para- 
graph 5-19* the general procedure of iteration method of design is given. 
The basic principle of the preliminary design of the column for any story 
after all the stories belov this story have been designed is given in this 
paragraph . 

The equivalent dynamic 
of a multi-story frame 
building is given in figure 
5-33- Any three consecutive 
floors, g - 1, g, and g + 1, 
are shown in figure 5.36. 

The equations of motion of 
these three floor masses are: 




a 2 * 



m 



'8 



g + 1 



dt 



g 



(5-79a) 



Figure 5.36. Three consecutive floors of a 
multi-story frame building 
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(5-79b) 



m 



g - 1 



- 1 



(5.79c) 



In the preliminary design of the columns for story g , it is assumed that 
all the stories above this story are absent. Neglecting R i n equa- 

tion (5.79*0; "the equation of motion for m becomes: 



m 



dt 



.(t) - 



15.80) 



> 



This is a typical equation of motion for a single-degree dynamic system. 
But the design method described in paragraph 5-12 is not applicable to the 
determination of R for one reason. In paragraph 5-12, the resistance 

O 

function depends on the absolute displacement of the mass supported by the 
columns. But in equation (5.80), although x is the absolute displace- 

o 

ment, the resistance function R depends on the relative displacement 

o 

between the floors g and g - 1 . Hence, in the preliminary design, it 
is necessary to consider the relative displacement between floors as well 
as the absolute displacement of each floor. 

The relative displacements, X , X^.-.X . ..X are defined by the 
following equations: 

X, = x. 



- 1 



X 



n 



x - x , 

n n - 1 



(5.81) 



Substituting equation (5.8l) into equations (5-79a)> (5-79^), & n4 (579c), 
the exact equations of motion in terms of relative displacements are ob- 
tained. By neglecting R in the preliminary design of columns for 

& J- 



story s , 'the equation of motion for mass, 
placement, X , is given by: 

ft 



where 



V 




in terms or relative dis- 



(5-82) 




X , and the resistance, 
&. 



Actiitl Eff*cli 
Ltti 



S3 

s 



: 




g 



Ic) Etlttlht 
Loon 



and f wv(t) is defined as the effective load for story g . The relative 
displacement, X , depends on this effective load which is a completely 

O 

knovn function since f n (t) and f (t) are the given loads and the 

g ~ ! g 

resistance function, R , is known from the preliminary design of the 

g - J- 

columns for story g - 1 . Once fy(t) is known, the relative displace- 
ment, X^, , and the resistance, R_ , at any time can be evaluated. 

A typical shape of the effec- 
tive load, ^(t) , is shown in 
figure 5-37a. The corresponding rela- 
tive displacement and resistance func- 
tion are given in figure 5.37b and 
figure 5-37c respectively. For pre- 
liminary design purposes, the actual 
resistance function shown in figure 
5-37c is replaced by an approximate 
resistance function shown in figure 
5-37d. The method of evaluating the 
parameters of the approximate resist- 
ance function, namely R . t 1 , t , 

g g ge 
and R is discussed in the follow- 

O^ 

ing paragraph. 

Linear Effective Load. 



(b) AttotWt 

Oiiplocmnl 



(c) Actual 
fttiiitanc 
Tune* ion 



Rtilitenct 

Fu " eliOB 



J 7 ^ 



b. 



In 





Figure 5.37. Relative displacement and 
resistance functions for a linear effec- 
tive load, no reverse yielding 



the preliminary design of two, three, 
or four story frame buildings, the 
initial portion of the effective load 
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for any story above the first can be approximated, with sufficient accuracy 
by a linear effective load as shown in figure 5-3?a. The initial portion 
of the actual x effective load is generally always negative because the initit 
notion of the lover floor masses is greater than the upper floor masses, 
jbis situation is reversed as the deflections increase. The general expres- 
sion of the linear effective load is: 



V*' 



B 



(5.83) 



Kbte that B ~ is a negative quantity for reasons mentioned above. Approx 



and T 



gX' 



nate formulas for the evaluation of the parameters, B 
given in paragraph 5-23- 

Substituting equation (5-83) into equation (5-82), the following 
equation is obtained: 




dt 



The resistance function, R , in equation (5-82) is replaced by kX in 

o o S 

equation (5-84). 

From equation (5.84), the time, t' , when X is zero, can be eval- 

o 

uated. The average value of R , that is R , in the time interval from 

/ 
to t ? and the relative velocity dX /dt at t' also can be evaluated 

g g 6 

These quantities, evaluated on the assumption that the columns in story g 

remain elastic when X is negative, are plotted in non-dimensional form 

O 

in figure 5-38- 

In actual design, a trial column section is selected from which the 
spring constant, k , the maximum elastic relative deflection, X ge , and 
the ultimate resistance, R , can be computed. The value of 
determined where T' is the period of a fictitious dynamic system ex- 



is 



pressed by equation (5-8 1 >-), and given by equation (5-85). 




(5-85) 



Using 



/T 1 in figure 5.38. the non-dimensional quantities: 
X gn 
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can be determined. From 
these non-dimensional quan- 
tities, the values of t 1 , 
_ . => 
R , and dX /dt can be com- 

o o 

puted. The value of t aa 
ge 

indicated in figure 537c can- 
be determined by assuming 
that the average relative ve- 
locity in the time interval 
from t' to t is 1.5 

CT fTP 
O O 

dX /dt , and this is given by: 
t = t 1 + ^r- 
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(5-86) 



, 

fi, linear effective load 
dt 



dt 



Once 



R 



and R 



^ , l^ . a.iju. 11 ou.c determined, the approximate shape of the 

D 6 O K^ 

resistance function at any time is known and this is shovn in figure 5-37d. 
Substituting R into equation (5-80), the maximum absolute displacement, 



gm 



g 
can be obtained by direct integration. Comparison of the computed 



- with the design specification determines the suitability of 



value of x_ 

the selected trial column section. 

In actual application, it is rather inconvenient to carry out the in- 
tegration process. In the following paragraph, a method is given for de- 
termining the suitability of the trial section based on use of the design 
charts presented in paragraph 5-10 for single-degree-of -freedom systems. 

oe 
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c. Equivalent System for Use with Design Charts. The design charts 
O f paragraph 5-10 vere constructed for the design of single-degree-of- 
freedom dynamic systems. These charts can be used to simplify the design 
of multi-story frame buildings by introducing an equivalent single-degree- 
O f -freedom system for the g-th story. 

In the previous paragraph the approximate resistance, R , at any 
time for the g-th story was developed. Suppose the external load, f (t), 
is a triangular load given by: 



(5.87) 



The load and resistance functions acting on m are shown in figure 5-39a. 

o 

The resistance function does not resemble the resistance function for a 
single-degree-of -freedom system with zero initial displacement and velocity. 
Therefore it must tie replaced by the resistance function, |E (t)l of an 
equivalent single-degree dynamic system shown in figure 5-39b. This is ac- 
complished by revising the external load in a manner consistent vith the 
change in the resistance function. For the equivalent system, the maximum 
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(a) Actual System 



(b) Equivalent Single-Degree 
System 



Figure 5.39. Equivalent single-degree dynamic system for the g-th story of a frame building 
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elastic deflection occurs is equal to t . The relationship between 

^ 



resistance is chosen to be equal to R and the time when the maximum 

elastic deflection occurs is equal 

r T 

| R g(t)L and R g (t) is given by: 

[ R g (t) ] eq =R g (t) - R> +R " (5-88) 



or 



where R' and R" are the curves of the approximate and equivalent re- 
sistance functions respectively for times less than t . R" is assumed 
to be a parabolic curve. 

Substituting equation (5-88) into equation (5.80), the following 
equation is obtained: 
2 

m &~2 " f g(*) - V*) = [ f g( t) ' R ' +R "] - [ R g (t) ]eq (5-89) 

Equation (5.89) indicates that the movement of the floor mass, m , can be 
represented by the movement of the mass of an equivalent single-degree dy- 
namic system. The external load on the equivalent system is: 

[f g (t)-R- + R"j 

instead of f g (t) . The work done on the equivalent system is: 

t 
TJ dt = / R " IT dt + fff (t) - R 2 ] ~ dt (I 

JQ J L g J dt 

The maximum allowable strain energy of the equivalent system is: 



dt gi 

where x ge is the absolute displacement of m at time, t 

t-> e 

90 
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The first term in equation (5.90) is identical with the first term in 
equation (5- 91)- Hence > in comparing the maximum work done and allowable 
energy absorption of the equivalent system, only the second terms in equa- 
tions (5.90) and (5-91) need be considered. The method described in para- 
graph 5-12* can be used for the evaluation of the maximum work done. The 
absolute maximum work done by the load Ff g (t) - Ril ^^ ^ flrst eval- 
uated. Next, figure 5-2? can be used to obtained the W0 rk done ratio, and 
the actual work done and allowable energy absorption determine the suit- 
ability of the selected trial section. 

It is seen from the above discussion that, by introducing an equiva- 
lent single-degree dynamic system, the design charts given in paragraph 
5-10 for single-degree dynamic systems can also be used for the design of 
multi-story frame buildings. The methods of evaluating the fictitious maxi- 
mum work done, work done ratio, and allowable energy absorption fo r tne 
equivalent system of the g-th story are described in the following three 
paragraphs . 



Ficitious Maximum Work Done, W . The fictitious maximum work done 



by the load 




(5.92) 



For a given value of T../T* > the values of R and t* can "be obtained 

gA gn g g 

from figure 5 '38- Within the range of T V /T' given by: 

gA gn 

0.15 < T v/T J < 1.0 
gX' gn 



the value of R ff ti varies from 

D G 



0.6 



to '0.85 



at ion may 
For the purpose of preliminary design, the following empirical eC & 

be used: 

H.,,0.7^ ^ (5 ' 93> 
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Substituting the value of R t obtained from equation (5-93) into equa- 
tion (5-92), the fictitious maximum work done, W^ , can be determined. 
Work Done Ratio, C The work done ratio, 

C, 



C W = W n/ W p as shown in 



figure 5.27 depends on C R and 
g-th story, C 



For the equivalent system of the 
can be approximated from the relation: 



R 



gm 



'6R 




(5.9"0 



but the ratio C 



gT 



has to be evaluated before C TT can be de- 

gw 



termined from figure 5-27- 

A single-degree dynamic system involves three non-dimensional param- 
eters. For the equivalent system of the g-th story, the parameters are 



chosen to be C 



v/v 



and 



C g T, 



where t 



ge 



is the time when the 



1.6 



1.4 




maximum elastic deflection, 

X , occurs. The relation- 
ge 

ship between these parameters 
is given in figure 5.^-Q for 
the case of a triangular load. 
This figure can be used for 
the evaluation of any param- 
eter when the other two pa- 
rameters are known. For the 
equivalent system of a g-th 
story, C _ is given by 

equation (5-9*0, t is 

ge 

given by equation (5.86), and 
T is the given load dura- 
tion. Hence, C _ and 

The C _ 



gR 
are known . 



Figure 5.40. t /T curves for elasto-plastic systems, 



f\s 144 CbO JI-" ClLl 
riangular load 



t ge /T g 

ratio can be obtained from 
figure 5-^0. The value thus 
obtained is defined as the 
equivalent value and denoted 



Using C gR &nd ( C g T ) 



eq 
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ratio, o can 



e eq , can 

tb en be determined from figure 5-2? from which the maximum work done, W 

determined. Prom figure 5-29, the time ratio, t ' 



can he 

and the time, 



when 



occurs can be determined. 



T is obtaind 



Maximum Allowable Energy Absorption. The maximum allowable energy 
absorption as represented by the second term in equation (5.91) is- 

/ v 



E 



R 



" x ge j 



(5-95) 



The value, 
given approximately by: 



x ge ' obtained fr m equation (5.80) by integration, is 




B t 




(5.96) 



tion 
C 



Comparison of Maximum Work Done and Maximum Allowable Energy Absorp- 

is given by equation (5-92). The work done ratio, 



The value -of 



, can be evaluated by the method described in the preceding paragraph. 
gw 

The maximum allowable strain energy, E , is determined from equation 

O 

(5- 95)- Comparison of the maximum work done and maximum allowable energy 
absorption determines the suitability of the selected trial section. 

Prom the above discussion, it is seen that, by introducing the equiva- 

lent system for the g-th story, and by using an equivalent value, /C \ , 

\ g j eq 
the design charts, figures 5-27, 5-28, and 5-29, which are prepared for 

single-degree dynamic systems, can be used for the preliminary design of 
the columns for a multi-story frame building. The detail procedure for de- 
sign is given in paragraph 5-23- 

5-22 ADJUSTMENT OF PRELIMINARY COLUMN RESISTANCE. In the design procedure 
outlined in paragraph 5-19, the preliminary design for each story is made 
assuming all higher stories are absent. That is while making a preliminary 
design of the g-th story, the presence of stories g + 1, g + 2, etc. is 
ignored. Actually, there is a force, R , acting on mass, m , from 

+ J. o 

the columns of story g + 1 (see figure 5.36). The effect of this force 
may be taken into account by an adjustment of the preliminary resistance of 
the columns for the g-th story. 
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The equation of motion for m , when both R and R 

B o e 

sidered, is given by equation (5-97) 
,2 



m 



g 



5-23 
+ i are con- 



- R 



(5-97) 



According to equation (5-97), the resistance function, R 



g + 



can be 



considered as an additional load on 
and 



m 



8 



The approximate shapes of R 



R as functions of time are shown in figure 5-41. When R 

Q o > -L 



g + 
i s 



neglected in the preliminary design of R , the maximum displacement, 

D 



R, 



f g + i^ 


R (g-H)m 


f (g+l)e 


1 



AR, 



gm 



Figure 5.41. Approximate resistance functions 



occurs at t . The problem now is 

to determine by what amount R 

gm 

should be increased so that when the 
additional load, R is acting 

on m , the maximum allowable dis- 

o 

.placement of m is still equal to 



gm 



The additional load, R 



o + 1 ' 

changes both the velocity and dis- 



placement of 
pose 



m 



at t 



gm 



Sup- 



R is changed by an amount 
equal to AR as shown in figure 

5.41; the magnitude of AR should 

gm 

be of such value as to cancel the ef- 



feet of R . 

g + 1 

equation (5-98). 



The modified value of R is given approximately by 



where 



Modified R 



gm 



gm ~ gm 




+ 1 



R gm = value obtain e<l in the preliminary design 
X - - V (V - V) - Vg 



g 



g 



R, 



+ t 
+ 1 g + 1 



(5.98) 



5 ' 23 COMPLETE DESIGN PROCEDURE FOR MULTI-STORY FRAME BUILDINGS. A complete 
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design procedure for the design of multi-story frame buildings is given in 
this paragraph. 

The procedure is set up on the assumptions that (l) the external 
loads can be approximated by triangular load curves, and (2) the actual ef- 
fective loads can be approximated by linear effective loads. The procedure 
is applicable only to the design of two, three, or four story buildings. 
The extension of this procedure to the design of higher buildings is dis- 
cussed in paragraph 5-24. Numerical examples illustrating the application 
of this procedure to the- design of two and three story buildings are given 
in EM 1110-3^5-^18- 

a. Preliminary Design of the First Story Columns. In the prelimi- 
nary design, the first floor is considered as a single-degree dynamic 
system with mass, m 1 , subjected to external load, f,(t) , and all the 
other floors are assumed to be absent. The energy method of design given 
in paragraph 5-12a is used. The maximum elastic deflection, X, , the 
plastic resistance, R , and the spring constant, k.^ , of the selected 
column section are determined. The time, t.^ , is determined from figure 
5.29 and the time, t. , is determined from figure 5>kQ. 

b. Preliminary Design of Columns Above the First Story. The proce- 
dure given in this section is applicable to the design of the columns of 
any story above the first story. Minor differences in the detail computa- 
tion for different stories are noted whenever it is necessary. 

Step 1. Evaluate the parameters, B and T , of the linear ef- 
fective load according to the following approximate formulas: 
For the second story: 



B 2X - ( m 2/ m l) B l - 



T, 




(5-99) 
(5.100) 



For any story above the second story: 

B =/m/m ,\ B i ~ B *j + R c l 
gX ^g'g-ljg- 1 8 B - " 



(5-101) 
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T g *( B - l)e B g - 1 g* 

where B I , 3 , B are the peak values of the external loads, 
f 2 (t), and f (t); and, T ] _, Tg, T g are the durations of the external 
loads. In determining B.^ Bg, B , T g , and T g , the given external 
loads should be approximated by triangular load curves. Depending on which 
story is under consideration, the approximation should be close in the fol- 
lowing time interval: 

for the second story, < t < 2t^ g , 

for the third story, 0-< t < 2tg e , and 

for the fourth story, < t < 2t. . 

In the design of buildings that are higher than two stories the value of 
B OT may be zero or negative in which case the procedure of paragraph 5-2^ 
should be used. 

Step 2. Compute R t' according to equation (5-93)- 

o o 

V T " 

R t' = 0.75 & 




Step 3 Compute the fictitious maximum work done, W , according 

. . ep * 

to equation (5-92). 




Before doing this, the given load, f (t) , is approximated by a re- 

D 

vised triangular load curve. Depending on which -story is under considera- 
tion, the approximation should be close in the following time interval: 
for the second story, < t < 1.3t , 
for the third story, < t < 1.3t , 
for the fourth story, < t < 1-St^ , and generally 

for the g-th story, < t < 1.3t/ N 
. (& ~ l)m 

The values of B g and T g obtained for the revised approximate tri- 
angular load will in general differ from those obtained in step 1 because 
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they simulate the given load over a different interval of time. 

Step U. Find R^ . As a first approximation, the following equa- 
tions may be used. 

for any intermediate story. R = W / x te. ->^\ 

S 111 gP gm O'-UoJ 

for the top story, R l.w /x 

gm J SP em. 

vhere x is the maximum allowable absolute displacement for the g-th 
story and is given by equations similar to equation (5.78), 

x lm = X lm 

v Y j. Y Y 
X 2m * A lm + A 2m X le 

V =s Y 4-Y J.Y Y V 

X 3m A lm + A 2m + X 3m " X le " X 2e 

Step 5- Estimate the maximum elastic deflection, X . If the 

heights of the first and g-th stories are equal, X may be assumed equal 

66 

to X, 



le 



Determine spring constant, k , from the relation, 

D 

R 



g gm' ge 
Step 6. Determine T 1 according to equation (5.85), 



(5-105) 



T' 2n 
gn 



/ m 

/p 7 8 

V t 1 + ("A - 1 



Step 7. Determine T _yA" where T is given by equation (5.100), 
or equation (5-102). 

Step 8. Using T /T' in figure 5 -38) find 



m 



m , / B 
g - I/ 8 



Then compute t 1 , R , and dX /dt . 

Q ^ 
Step 9. Compute t according to equation (5.86), 



ge 



ge 



Step 10. Compute t /T and 



from equation (5- 
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R /B 

c s 



R 



_ 

VB 

Step 11. Using t ge /T g , and C^ , determine (C gT ) eqL from figure 

5.UO. 

Step 12. Using C^ and (c^)^ , determine C gW and t^/Tg from 
figures 5.27 and 5-29 respectively. Then compute W^ and t^ from the 
relations, 

V - W (5 ' 106) 

(5 ' 107) 



The value of t is used to check whether the triangular load approxima- 

gm 
tion in step 3 is acceptable. That is, if t is much different from 

1 3t, \ as assumed in step 3, a revision of the triangular load ap- 

(g - Ijm 
proxiraation may be necessary. 

Step 13. Compute x according to equation 




Step Ik. Determine the maximum allowable strain energy, E , ac- 

O 

cording to equation (5-95). 

g gm \ gm ge] 
Step 15. Comparison of E and W determines the suitability of 

R . If the agreement between E and W is not within 10$, the pro- 
cedure from step 3 to step ik is repeated by using a new value of R 
given by, 

E + W 



R 



gm 2 x 



- ^ 
x - x \ 
\ gm gej 




However, if E and W agree within 10$, the value of R obtained 
' g gin ' ' gm 

from equation (5-103) is approximately the maximum resistance required. 

c. Adjustment of Preliminary Column Resistance R . The effect of 

gm 

the g + 1 story columns on the movement of floor mass, m , is taken into 

' g ' 

98 
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account by modifying the value of 
given by equation (5-98), as 
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The modified value of R I B 

gm 




where 



Modified R = R 



R = value obtained in the preliminary design 



= 






V - *ge) 



L g + 1 , E 






(g -t- 



' R g + l*i + 1 



Starti-na from the top story, the plastic resistances of different stories 
should be modified one after the other. In the modification of R by 
equation (5-98), the value of R^ g+ ^ should be the modified value 
rather than the value obtained in the preliminary design. In the computa- 
tions of I-, for the bottom story, the value of R. is zero and the value 

of t is equal to 2t, /3 . 
ge ie 

d. Selection of Final Column Sections. The maximum horizontal re- 
sistance of the selected columns for the g-th story, including the effect 
of vertical loads acting on the columns as outlined below, should be close 
to the modified value of R obtained in the previous paragraph. In de- 
termining the effect of vertical load, its average value in the time 

interval, 

t < t < t 
ge gm 

should be used. The expression of the ultimate resistance of the selected 
columns is given by: 



p _ 

gm 



2h 



(5-109) 



where 



M^ = plastic resisting moment of the selected section after the ef- 

feet of the vertical load is taken into account. 
ZP . the time average of the total vertical load in the time interval, 

v t < t < t 

ge gm 

n = number of columns in the g-th story. 
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h = height of the g-th story. 

After the column section has been selected, the spring constant, k , 
and the limiting elastic deflection, X . are determined. In the evalua- 

o e 

tion of these parameters, the effect of girder flexibility may be approxi-. 
mately accounted for (see EM 1110-3^5-^17 and EM 1110-3^5-^18). 

e. Numerical Analysis of the Entire Building. A numerical analysis 
for the entire building should be carried out to determine the suitability 
of the selected columns. If the result of the numerical analysis indicates 
that the maximum absolute displacement, x , is greatly different from 
the design specification, all the columns in and below the g-th story 
should be revised. The revised value 
and error method or from the relation, 



should be revised. The revised value of R can be determined by a trial 



p 

= 



" R (g + pmfen " X ge\ 



R (g + l)m( X gm 



x -X 

gm ge 



(5-110) 



where R 1 and R5 .. \ are the average values of the maximum resistance 
gm (g + IJm 

used in the numerical analysis over the time interval, t < t < t ; and 

x' and x' are the displacements obtained from the numerical analysis. 
gm ge 

R/ ., v is the revised maximum resistance of story g + 1 : x is the 
(g + l;m gm 

prescribed maximum allowable absolute displacement; and f It' j is the 
magnitude of the external load, f (t) , at the time when x 1 occurs. 

> &IH 

New column sections should be selected according to the revised 
values of the plastic resistances. The spring constants and the maximum 
elastic deflection should then be determined. The numerical analysis 
should be repeated using the revised column sections. 

5-2U PROCEDURE FOR BUILDINGS OF MORE THAN TWO STORIES. The design pro- 
cedure given in paragraph 5-23 is for the case of a linear effective load 
which is applicable only to buildings two stories or less in height. In 
the design of higher buildings (higher than two stories) the value of B Y 

-A 

computed from equation (5-99) m y be zero or negative. The effective load, 
' and the relative displacement, X g , are then positive immediately 
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B computed from equation (5-101) is zero or negative, the value of t 
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after time, t = . Thus when B^ computed from equation (5.99) is zero 
or negative, the time, tg e , can be taken to be 5t le / 3 . Similarly, if 

b 3e- 

can be taken to be I t ge + -~l . For buildings of more than two stories, 
if B^ is close to zero or negative, R g should be considered equal to 
zero and the procedures from step 5 to step 9 in paragraph 5-23b should be 
omitted in the preliminary design of the second or third story columns. 
5-25' ACCURACY OF FIOOR-BY-FIflOR DESIGN METHOD. A number of assumptions 
have been used in this design method. First is the replacement of actual 
effective load by a linear load in which the parameters, B and T 
are determined by approximate formulas. Second, approximate formulas are 
used in the determination of t and C . Third, in the evaluation of 

O 1 - &T\ 

strain energies and work done, the portion of R"(t) shown in figure 5-39b 
is neglected. And fourth, in the determination of the plastic resistance, 
R , of the chosen column section from equation (5.109), the time average 
of the total vertical load on the columns is used. 

There is error involved in each assumption. The accumulative error 
for an actual design problem may be quite large. However, based on a 
number of examples which have been worked covering a wide range of build- 
ings under various loadings, the maximum resistance of the columns for any 
story obtained by this floor-by-floor iteration method has been found to 
be generally within 20$ of the desired value, for cases in which the pre- 
scribed maximum allowable relative displacement was at least three times 
larger than the maximum relative elastic deflection. 

For multi-story frame buildings, a small change in the plastic re- 
sistance generally results in a large change in the relative displacement. 
Although the result of the numerical analysis may indicate a large error 
in the relative displacement, this usually represents small error in the 
plastic resistance. Since the plastic resistance cannot be evaluated 
exactly, accuracy greater than 20$ in the design of columns on the basis of 
relative displacement is not warranted. An error of 20$ in the relative 
displacement may correspond to an error of only 5$ in the plastic 
resistance. 
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DESIGN OF MULTI-STORY BUILDINGS FOR ELASTIC DEFORMATION 

5-26 INTRODUCTION. In the design of a multi-story building for elastic 
deformation the maximum allowable lateral deflection of the columns, under 
dynamic load, should not exceed that which causes yielding. In the equiva- 
lent system for frame buildings as shown in figure 5-33> the moment in the 
columns of the g-th story is proportional to the relative displacement X . 

O 

Direct stresses are developed in the columns due to over-all bending but 
these may be ignored in buildings of three or four stories. In frame 
buildings, where the direct stress due to over-all bending is appreciable, 
the preliminary design could still be made without including the effect of 
the direct stress. This would be considered, however, in selecting the 
column sections. Hence, the elastic design of frame buildings for elastic 
deformation can be based on the requirement that the maximum allowable 
relative displacement for any story must be smaller than the relative 
elastic deflection at the yield point stress. 

However, in the equivalent system of a shear wall building as shown 
in figure 5.32, the stresses in the shear wall at any story are a function 
of the over-all bending and shear in the wall and hence depend on the rela- 
tive displacements of all stories rather than just the relative displacement 
of adjacent stories. The use of absolute displacements is desirable in the 
design of shear walls because it is the simplest method to express the rela- 
tive displacement between all floors in a shear wall building, whether the 
floors are adjacent or non-adjacent. 

The characteristics of a multi-story building designed for plastic 
deformation discussed in paragraph 5-20 are also generally true when de- 
signed for elastic deformation. These characteristics are: (l) the rela- 
tive displacements in the upper stories are negative shortly after the 
initiation of the load; and (2) the elastic rebound of floor masses must 
be considered in the evaluation of the relative displacements. In design 
for plaatic deformation, the resistance, R , is assumed to be a constant 

O 

after the elastic limit is reached. The design of the g-th story columns 
is reduced to the determination of the shape R as shown in figure 5-37^- 

O 

Moreover, the elastic rebound of floor masses is considerably smaller than 
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tne allowable relative displacement. Hence, the use of approximate formula 
in estinating the magnitude of elastic rebound is permissible. But in an 
elastic system, the resistance, R g , always varies with time. The elastic 
rebound is about of the same order of magnitude as the relative displace- 
aent, and approximate evaluation of the elastic rebound is not permissible. 
Because of these differences, the elastic design of multi-story buildings 
is more involved than the design for plastic deformation. In most cases, 
a trial and error design procedure may have to be used. The elastic de- 
sign of multi-story frame and shear wall buildings is briefly discussed 
in paragraphs 5-27 and 5-28 respectively. The design of shear wail build- 
ings is more fully discussed in EM 1110-3^5-^19. 

5.27 DESIGN 0? MULTI-STOKY FRAME BUILDINGS FOR ELASTIC DEFORMATION. One 
rjossible approach -to the elastic design of multi-story buildings is the 
preparation of non-dimensional graphs. The feasibility of this method is 
illustrated by the characteristics of a two-story frame building. The 
equivalent dynamic system of the 
building is shown in figure 
5^2. The external loads, f -^ 

and f (t) are assumed to be 

2 
triangular with load durations 

T and T respectively. The 

relative displacements, X ls B| ( I - t /T, ) = f, ( t ) 

and X- when the peak external 
2s 

loads are applied statically, 
are given by 




/77T7 



*2s - 2k, 



Figure 5.42. Equivalent dynamic system for 
two-story building 



B 



B 



X 



Is 



2k, 



(5-112) 



Suppose 



and X denote the maxinrum relative displacements under dy- 



.^ 



^ 
namic load^ the dynamic load factors for the second and first stories are 

defined by the following relations: 

fn T v - (5-H3) 

(D.L.F. 
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5-2? 



These dynamic load factors depend on the non-dimensional parameters 

T 2 /T 1' B 2/ B 1> and V T ln ' Wbere Tln ^ 6C1Ual t0 2 * *V k l ' "*** 6en- 
eral shape of the dynamic load factor curves is shown in figure 5-^3 for 

specific values of kg/k.^ Tg/T^ 
and B 2 /B, . Before these curves 
can be used for practical design, a 
series of similar graphs must be pre- 
pared for different combinations of ' 

m^/m,, B 2/ B i and T 2/ T 1 to cover 
the probable range of practical 
buildings. For practical cases, it 
is believed the range of nu/m., 

B 2/ B 1> T 2/ T 1' and ^^1 W uld vary 
within narrow limits. 




k z /k| incrtailng 



I? 




Incrtailng 



T, /T., 



Fig 5.43. Dynamic load factors for too-story 

building for elastic design, nvj/m^ = constant, 

r 9 /T = constant, B^l ~ constant 



It is even 

possible that they might be assumed 
constant for certain types of con- 
struction. Suppose four different 
values are used for each non- 
dimensional parameter, then the 
total number of graphs needed is 6k. 
By similar considerations for a 

three-story building, the total number of graphs required is 2^01. Since 
such a large number of graphs are needed, it appears impractical in general 
to prepare non-dimensional graphs for use in the design of multi-story dy- 
namic systems. 

For practical design, a trial and error procedure may have to be 
used. A set of columns for the entire building is selected. The suit- 
ability of the selected columns is determined by the result of a numerical 
analysis. In the preliminary selection of columns for the g-th story, an 
equivalent single-degree dynamic system, shown in figure 5-^, may be used. 
In the equivalent system, all the floor masses and external loads above 
the g-th story are assumed to be concentrated at the g-th floor. The 
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...m 



Figure 5.44. Equivalent single-degree system for the g-th story 

design method on the basis of dynamic load factor given in paragraph 5-12c 
is used for the determination of the columns for the g-th story. The 
system shown in figure 5-^ is used for the preliminary design of columns 
in the first, second, up to (n - l)th story. For the top story, the ex- 
ternal load is assumed to he 2f (t) instead of f (t) . 

Because of higher modes of oscillation in multi-story buildings, a 
large number of maximum and minimum values are presented when the relative 
displacements are plotted versus time. Hence, in the numerical analysis, 
the integration process should be carried out far enough to insure that 
the absolute maximum relative displacements for all stories are obtained. 
5-28 DESIGN OF MULTI-STORY SHEAR WALL BUILDINGS FOR ELASTIC DEFORMATION. 
In the elastic design of shear wall buildings, the controlling factor may 
be shear and/or bending stresses. For the evaluation of both stresses, the 
building is replaced by a multi-degree dynamic system discussed in para- 
graph 5-l6a. The method of evaluating the parameters of the dynamic system, 
the detailed preliminary design procedure, a discussion of design tech- 
niques for both elastic and plastic action, and numerical examples of 
shear wall building design are given in EM 
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APPENDIX A 



EFFECT OF VARIATION OF PARAMETERS ON THE RESPONSE 
OF SINGLE-DEGREE DYNAMIC SYSTEMS 



A-01 INTRODUCTION. In paragraph 5-10, the response of single-degree dy- 
namic systems subjected to simplified loads is presented in non-dimensional 
graphs. For & given problem, when the peak load B , the load duration T 
the mass m , the spring constant k and the plastic resistance R are 
all known, the maximum work done, W ffl , the maximum displacement, x , and 
the time, t m , when x ffl occurs can be determined from these non- 
dimensional graphs. For the convenience of discussion, the first five var- 
iables (B, T, m, k, and R m ) are considered as independent' variables in 
this appendix. 

For design applications, in addition to these non-dimensional graphs, 
it would be helpful to know how the response of a dynamic system varies 
with the shape, the peak value, and the duration of the load pulse, and how 
the response is affected by the mass, the spring constant, and the plastic 
resistance. In other words, it is of interest to know the effects of small 
increments of the independent variables on the response. Two practical de- 
sign problems, in which this is of importance, are the following: 

The first problem is concerned with uncertainties. Since all the in- 
dependent variables are subject to uncertainties in their determination, it 
would be of interest to know what is the effect of these uncertainties on 
the response and what is the necessary accuracy in the determination of the 
variables such that the response of a dynamics system is within a given 
tolerable limit. 

The second problem occurs when it is intended to answer the following 
two questions which often arise a.t the completion of a trial design compu- 
tation. These questions are: (l) What is the permissible difference be- 
tween the specified maximum displacement and the corresponding computed 
value of a given design; and (2) if a trial design is judged to be .unsatis- 
factory, by what amounts should the variables be changed in the next trial 
design. For numerical examples, suppose for a given system a 10$ change 
in the mass results in 1# change in the maximum displacement, then an 
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accurate determination of the mass would not be necessary. On the other 
hand; if a 10$ change in the plastic resistance results in a 50$ change 

in x , a close agreement between the specified maximum displacement and 

m 

the corresponding computed value would not be warranted since the deter- 
mination of the plastic resistance may be subject to an uncertainty greater 
than 10$. 

The effects of small increments of the independent variables on the 
response are discussed in this appendix. These effects are expressed in 
terms of percentage-increment ratios which are non-dimensional quantities 
and are defined in paragraph A-02. In paragraph A-03, the technique of 
non-dimensional analysis is applied in deriving certain functional rela- 
tionships among the percentage-increment ratios. Because of these rela- 
tionships, the effect of small increments of all of the variables (B, T, 

m, k, and R ) can be condensed into two sets of curves. The percentage- 

m 

increment ratio for linearly elastic, completely plastic, and elasto- 
plastic systems subjected to both rectangular and triangular loads are dis- 
cussed in paragraphs A-OU, A-05, and A-06 respectively. The effect of load 
shape on the response is discussed in paragraph A-07- The expressions of 
various percentage-increment ratios for linearly elastic, completely 
plastic, and elasto-plastic systems subjected to rectangular and triangular, 
loads are summarized in paragraph A-08. 

It is to be noted that all of the material in this appendix is de- 
rived from the corresponding material in EM 1110-3^5-^15 through a process 

of differentiation. The percentage^ increment ratios are related to the 

oXj^ oXj,. 

partial derivatives, (for example, -^ , -^)j and the graphs in this ap- 
pendix are related to the slopes of the corresponding graphs in paragraph 
5-10. Although these two sets of graphs are derivable from each other, 
they are intended for different applications. Generally speaking, for 
routine design computation, it would be easier and the result more accu- 
rate when the graphs in paragraph 5-10 are used, while for those special 
applications discussed earlier in this paragraph, it would be better to 
use the graphs in this appendix. Beside these special applications, the 
graphs of percentage -increment ratios also illustrate many important prop- 
erties of structures under dynamic loads . For this reason, although , 

ion 
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quantitative data are presented in figures A.I through A. 8, the discussion 
In paragraphs A-(fc, A-05, and A-06 is more or less qualitative in nature. 
Jhis quantitative discussion should then serve as a rule of thumb for a 
designer to make proper judgement on the many practical problems which are 
not necessarily related to any routine design computation. 
A-02 SENSITIVITY AND HERCENTAGE-INCREMENT RATIO, a. Sensitivity. The 
maximum displacement, x m , generally depends on all the variables, B, 
T, m, k, and R m . Suppose for the time being only the load duration, 
T , is varying. Then x ffi is a function of T . If the load duration is 
changed from T to T + AT , the corresponding maximum displacement is 

changed from x to x + Ax . For small increments, Ax and AT are 
m m in in 

related by equation (A.l): 



o-x 



'm 



~5T 



AT 



(A.I) 



or 



xT 



Ax 

m 

AT 



(A. 2) 



where 



AT s the increment of the load duration. 
Ax = the corresponding increment of the maximum displacement. 

o"x 

-ip = a partial derivative. 

The quantity, S , being a partial derivative, is defined as the sensi- 

A JL 

tivity of the maximum displacement with respect to the load duration. Thus 
the sensitivity S _ is equal to the ratio of two increments Ax and AT 
and is a dimensional quantity: For example, when S XT is equal to 0.01 
inch/sec, this indicates that if the load duration is increased by 0.1 
seconds, the maximum displacement would be increased by 0.001 inches pro- 
vided the increments are small in comparison with their original values. 
Although the dimensional quantity, sensitivity, is very useful in 
numerical computation, it does not express clearly the physical significance 
of "sensitivity." From the numerical magnitude of the sensitivity, it is. 
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difficult to perceive whether the maximum displacement is "very sensitive" 
or "not sensitive" to any variation of the load duration. This is because 
the original values x and T to which Ax m and AT are to be added 
are not included in the definition of the sensitivity. Hence in this ap- 
pendix, instead of using increments and sensitivities, the non-dimensional 
quantities, percentage-increments and percentage -increment ratios are used. 
These non-dimensional quantities are defined in paragraph A-02b. 

b. Percentage- Increment Ratios. If the left hand side of equation 
(A.l) is divided by x m and the right hand side is divided by -^ T , the 
result is given by equation (A.3). 




= I xT (A.3) 

where . 

Ax 

= percentage-increment in x 



m 



AT 

- = percentage-increment in T 




and 

i -i 

XT x 

m 

= The ratio of the percentage change in x to the percentage 
change in T while all the remaining independent variables are 
held constant. 

The quantity I T , with two subscripts, is defined as the percentage- 
increment ratio. 

In the discussion so far, x and T are considered respectively 
as the dependent and independent variables. By a similar procedure, if the 
variables are the maximum work done and the peak load B , the percentage- 
increment ratio L^ would give the ratio of the percentage -increment of 
W ta to the corresponding percentage -increment of the peak load B . 

There are a number of percentage-increment ratios for a single-degree 
dynamic system subjected to simplified loads. These ratios are described 

in paragraph A-02c. 
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c ' Various Percentage-Increment Ratios for Single-Degree Dvnam^ 
Systems^ For a single-degree dynamic system subjected to simplified loads 
there are five independent variables, _B, T, m, k, and R . when any 
one of these variables is changed, the maximum displacement, in general, is 
affected. The total effect on x m when all the variables are changed can 
be obtained approximately by the method of superposition. Thus the total 
increment Ax , for the general case instead of being given by equation 
(A.l), is given by: 



m 



AB + 



m 



AR 



(A. 10 



m 



According to equation (A. 2), the different partial derivatives in equation 
(A.I*) can be defined as different sensitivities. 

Applying the technique used in deriving equation (A. 3), the general 
expression for the percentage -increment in the maximum displacement can be 
derived and this is given by equation (A. 5). 



Ax 



m 



m 



(A. 5) 



m 



xT T ' "xm m T x xk k T x xB B " ^cR IT 
By the same method of derivation, if W instead of x_ is considered as 



m 



m 



the independent variable, the general expression of the percentage- increment 
of the maximum work done is given by equation (A.6): \ 



AW 

-J5 _ T 
W m ~ J VTT' 



AR 



The quantities, I 



'XT' 



xm 



1. 
Stem 



xk' 



m 



T s^: j. T - CA\ fi ^ 

"Sffi B + -Sffi R (A '\ } 

\ 

and !.. in equation (A. 5) are tlae 



"xR 



percentage -increment ratios for the maximum displacement, while the quant i-\ 
ties, 1^, 1^, !^ f i^ and 'l^ in equation (A.6) are the 
percentage-increment. ratios for the maximum work done. 

The percentage-increment ratio gives quantitatively the effect of a 
given variable on the response of the dynamic system. For example if the 

numerical values of I m for different structures are' different, the 

xT 
larger the value of I is, the more the maximum displacement will be af- 

xT ' 
fected by any change in the load duration. Similarly, for a given structure 
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if I is numerically larger than 1^ , the maximum displacement will be 
affected more by the plastic resistance than by the mass, m . On the other 

hand, if I , is numerically small, say less than 0.2, the maximum dis- 
J xk 

placement is practically a constant when the spring constant k varies. 

Equations (A. 5) and (A. 6) are derived by partial differentiation and 
superposition. Hence, in the evaluation of the percentage change in x 
or W . for a given structure dynamic system, these simplified equations, 
in general, are applicable only when the increments in B, T, m, k, or 

R are not more than 10$. 
m 

In equations (A. 5) and (A. 6), there are ten percentage-increment 
ratios. These ratios are not all independent. The relationships among 
them are derived in paragraph A-03- 

A-03 RELATIONSHIPS AMONG THE PERCENTAGE -INCREMENT RATIOS, a. Increment 
Ratio for Non-Dimensional Parameters. The discussion in paragraph 5-10 
shows that only three non-dimensional parameters are needed to express the 
functional relationships of all the variables of an elasto -plastic system 

subjected to simplified loads. For example in figures 5-25 and 5-28, the 

x R 

three parameters are, , 7^- and -^- , which are denoted by C , C , 

x e i n a xT 

and CT, respectively. That is 

K 

C =^! (A. 7 ) 

xx \ i / 

e 



In terms of these parameters, the dependent variable of the dynamic system 
is C^ instead of x m , and the independent variables are C^ and C_ 
instead of B, T, m, k and R m . For small increments in C and C R , 
the corresponding increment C is given by equation (A. 10). 

A 



112 



A ' 3a 



EM 1110-31*5-415 
15 Mar 57 



vhere M^, AC T and AC R are small increments, j and 
tial derivatives. T 

Applying the same technique used in the derivation of equation (A. 3) 
and (A. 5), the percentage-increment of G X , can be derived and is given- by 
equation (A.ll) 



where 



~ + J 2 T (A - U) 




=- = percentage -increment in C 
C x 

C T 

s- = percentage-increment in C 



= percentage -increment in 



R 



T 
the percentage-increment ratio of C with respect to CL 

X 



T 
2 = C x 



= the percentage- increment ratio of C with respect to CL . 

x R 

There are only two percentage-increment ratios in equation (A.ll). 
It can be concluded that of all the percentage-increment ratios of an 
elasto-plastic system, only two ratios are independent. Hence the ten 
percentage-increment ratios in equations (A. 5) and (A. 6) can be expressed 
in terms of 1^ and 1^ . 

This conclusion can also be deducted by physical reasoning. C is 

X. 



a function of C and C m . If for a given case C is a constant, C 

.n J. r\ 

would remain a constant provided that C_ is unchanged although T , m 
and/or k may be varying. This indicates that I _, I , and I . can 

XJ. * 



XIU 



be expressed in terms of I, . Similarly, it can be concluded that 



, . 



an <l 



can be expressed in terms of I 2 . The equations expressing 
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those relationships ore given in paragraph A-03b. 

b. Percentage-Increment Ratios in Terms of I- L and I g Substi- 

R 

tuting -=2 for x in equation (A.?), the following equation is obtained: 
K. c 

x k 

C,i (A.*) 






Tfeking logarithms of both sides of this equation, then differentiating, and 
then replacing the differential notation "d" by the increment notation 
"A," equation (A. 13) is obtained. 

5i ^m + Ak f\n (A } 

C x~ X m k R m 

Similarly, from equations (A. 8) and (A. 9), the following two equations are 
obtained respectively: 

AC 

T _ AT _! A^ 2: ^JB (A "iM 

~C~ = ~T~ + 2~~2m 

AC D AR AB 

1 _ A (A 15) 

C ~ R B ^ P ' 
R m 

substituting equations (A. 13), (A. 1*0 and (A. 15) into equation (A. 11) and 
rearranging: 

Ax I/I\ R / ^ 

m \ / m 

Comparing equation (A.l6) with equation (A. 5), the following relationships 
among the percentage -increment ratios are obtained: 

(A. 17) 



So far in the derivation, x is considered as the dependent variable. By 

' 

emen 
111:- 



the same method, the percentage -increment ratio when W is the dependent 
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variable can be derived. For an elasto-plastic system, referring to equa- 
tion (5-38)> the expression for W ffl is given by 



(A.19) 



Taking logarithms of both sides, differentiating and rearranging: 



> 



m 



2C x -l 



m 



2C x -l 



'm 



L Ak 

-Ik. 



(A. 20) 



Substituting equation (A.l6) into (A. 20), the expression of -~ in terms 
of the percentage changes of the dependent variables is obtained: 



2C 



A m 
AT 

T 



Am 
m 




2C - 1 ^ T + ( 2 + 2C - 1 ^l "R~ 
x \ x / m 



Comparing equation (A. 21) with equation (A. 6), the relationships among the 
percentage-increment ratios < 
equations (A. 22) and (A. 23). 



percentage-increment ratios of W are obtained and these are given by 



2C 



1 - 2C 2C - 1 
x T x 

C -Sto " C 



1 - 2C 



2C 



2C - 1 
x T _ x 

"Sffi ~ 2C 



(A. 22) 
(A. 23) 



Equations (A. 20) through (A. 23) are derived on the basis that W ffi is 
given by equation (A. 19), which is true for an elasto-plastic system. It 
can be shown that these equations are also applicable for both linearly 
elastic and completely plastic systems provided that C^ is replaced by 
unity for the former while the limiting value with G X approaching in- 
finity is used for the latter. 

Equations (A.I?) and (A. 22) show that 1^ , 1^ , 1^ , 1^ t ^ 
and L^ can be expressed in terms of 1^ , while equations (A.18) and (A. 23) 



show that I 



and 



can be expressed in terms of 



XB XR nJJ "*> , . . 

Therefore, the graphical presentation of the percentage- increment ratios 
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is greatly simplified by these functional relationships. When any one of 
the- percentage-increment ratios in the first group is computed, all the re- 
maining I's are also known. The same is true for the second group of 
I's . Thus only two sets of curves are needed for the complete graphical 
presentation of all the I's of an elasto-plastic system subjected to 
simplified loads. In this appendix, one set of curves is for 1^ , which 
is also equal to I T . The other set is for I _ which is equal to 

(. X J.' Jtr\ 

1- + I. J . These curves are described in paragraph A-06. 

c. Other Percentage-Increment Ratios. The discussion so far covers 
the cases when the independent variables are B, T, m, k, and R ; and 
the dependent variables are x and W ffl . For special application, the 
percentage-increment ratio for variables other than those mentioned above 
may be needed. Because of the non-dimensional nature of the percentage- 
increment ratio, the expressions for additional percentage -increment ratios 
can be derived from those already presented by a method of superposition. 
A few examples are given for illustration. 

As a first example, suppose the effect of the limiting elastic dis- 

placement x on x is to be considered. Since 

m 

R 
m 

x e = T 



therefore 



Ax AR ., 

_ e __ m Ak 

x ' R k 
e m 



The percentage-increment ratio of x with respect to x is then given 

me 



by: 

T xe 

Where I xe is ratio of the percentage change in x to the corresponding 

percentage change in x . 

e 

As a second example, the load pulse H , which is the area under the 
load time curve as defined by equation (5-5), is considered as an inde- 
pendent variable. For rectangular or triangular load, the load pulse is 
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proportional to the product of B and T 
obtained: 

I = I _ + 
xH xB 
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Then equation (A. 26) is 



(A. 26) 



where 1^ is the percentage-increment ratio of the maximum displacement 

x with respect to the load pulse H . 

m B 

As a third example, the effect of = on the maximum displacement is 

considered. For a given triangular or rectangular load, this variable is 
the ratio of the two sides of" the triangle or rectangle. Hence in a very 
restricted sense, it represents the "shape" or "proportion" of the load 
curve with a given geometry. Let this variable be denoted by J , that is: 



(A.2?) 



The 



percentage- increment ratio I T is given by: 

xj 



T - T - T 
XJ XB XT 



(A. 28) 



Where I T is the ratio of the percentage change in x to the correspond- 

B m 

ing percentage change in 7= - 

A-04 PERCFJiTAGE-INCpMEM' RATIOS FOR LINEARLY EIASTIC SYSTEMS. The 
percentage -increment ratios for the maximum displacements of linearly 
elastic systems subjected to triangular or rectangular loads are given in 
this paragraph. In deriving the expressions for different percentage- 
increment ratios, the methods described in paragraphs 5-07b, A-02 and A-03 

T 
are used. For example, the case of a triangular load with ^- > 0-371 is 

considered. This corresponds to t^ < T . The expression for x m for 



m 



this case is given by equation (A. 29). 



B 

x = r 
m k 



(A. 29) 



The expressions for I n and 



xB 



spectively given by the following two equations: 



I derived from this equation are re- 

xT 



XB 



1.0 



(A-30) 
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xT 



tan 2C T 



JL _ 1 



. 31) 



Applying equation (A. 17), the corresponding expression for I and I 

"^ Ll A-Ix 

can be derived from equation (A.3l). 

Similar expressions for all the percentage-increment ratios for tri- 

m 

angular loads with ~ < 0-371 , and also for rectangular loads can be ob- 
tained. These expressions are all summarized in paragraph A-08a. 

The percentage- increment ratios for linearly elastic systems are 
plotted in figures A.I and A. 2. The first figure is for rectangular load 
while the second is for triangular load. From these two figures, the fol- 
lowing conclusions can be drawn: 

(1) For both types of load, I is a constant and equal to unity. 

XD 

The maximum displacement of a linearly elastic system is always directly 
proportional to the peak load B . 

(2) When C T is less than 0.2, I is approximately constant and 
equals unity for both types of loads. Thus for short duration load, the 
maximum displacement of a linearly elastic system is directly proportional 
to the load duration T . On the other hand, for a rectangular load, I ^ 
is zero when C T > 0.5 and for a triangular load, I is less than 0.1 

L JS. JL 

when C T > 3-0. Therefore the maximum displacement is practically inde- 
pendent of the load duration for long duration load. 

(3) For both types of loads, when C < 0.2, both I and I are 

X X1U XlC 

approximately equal to -0.5- Thus for rectangular or triangular loads of 
short duration the maximum displacement is inversely proportional to both 
the square root of the mass m and the square root of the spring constant 
k . On the other hand, for a rectangular load when C T > 0. 5, I equals 
zero and I . equals -1.0. These numerical values are also approximately 
correct for a triangular load when C > 3.0. Thus for long duration tri- 
angular 'or rectangular load, the maximum displacement of a linear elastic, 
system is- practically independent of the mass but is inversely proportional 
to the spring constant k . 
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Figure A.I. Percentage increment ratios for a linearly elastic system 
subjected to rectangular loads 



1.5 



o 

K 

c 
II 

E 



Trlonqulor Resistance -, 




-1.5 



A2. Percentage increment ratios for .linearly elastic system 
subjected to triangular loads 
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The numerical values of 1^, I XT , 1^ and 1^ axe always 
leaa than unity. Thus the reliability of the maximum displacement de- 
termined by a design computation is better than or as good as the uncer- 
tainties in the determination of any variable. 

A-05 FCTCENTAGE-INCREMENT RATIOS FOR COMPEETEU PIASTIC SYSTEMS . Applying 
the method described in paragraph 5-07c, the expression of x m for a com- 
pletely plastic system subjected to rectangular loads can be obtained and 
is given by equation (A. 32). 






The percentage -increment ratios 1^ and 1^ derived from equation 
(A. 32) are given by: 




(A.33) 



The corresponding expressions for 1^, 1^, 1^, 1^, 1^ and 
can be obtained from equations (A. 17), (A.18), (A. 22) and (A. 23). Similarly 
the corresponding expressions for triangular loads can be derived. These 
expressions are summarized in paragraph A-08b. 

The percentage-increment ratios for both x and W of a completely 
plastic system subjected to rectangular and triangular loads are given in 
figures A. 3 and A A respectively. Prom these tvo figures, the following 
conclusions can be drawn: 

(1) For both types of loads, I XT and I are equal to 2, and I 

and 1^ are equal to unity. Hence, for a completely plastic system sub- 
win 

Jected to simplified loads the maximum displacement x and the maximum 

m 

work done W m are: (l) proportional to the square of load duration T 
and (2) inversely proportional to the mass m . 

(2) The numerical values of 1^ or 3_ are always greater than 2 
and the numerical value of 1^ is always greater than unity. When 

Cp < 0.3, 1^ and 1.^ are approximately equal to 2 and -1 respectively. 
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,rttire A.3. Percentage increment ratios for a 
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Figure A.4. Percentage increment ratios for a 

completely plastic system subjected to 

triangular loads 



re- 



Hence, for such cases the maximum displacement is proportional to the 
square of the peak load, B , and inversely proportional to the plastic 

sistance, R . 
m 

(3) When C R > 0.5, for rectangular loads, I > 3 and I < -2; 
and for triangular loads, I > 4 and I < -3. Thus the maximum dis- 

nl 4- ^^ ^^ 

Placement is extremely sensitive to variations in B and R . For such 

m 
cases it would be meaningless, in a trial and error design procedure, to 

seek an agreement better than +30$ between the maximum displacement obtained 
by a design computation and that of the design specification. 
A ~6 PERCENTAGE- INCREMENT RATIOS FOR ELASTO-FLASTIC SYSTOMS. The expres- 
sions for the percentage -increment ratio of an elasto-plastic system sub- 
jected to simplified loads can be 'derived by the methods described in para- 
graphs 5-0?d, A-02 and A-03- An example is used to illustrate the procedure 
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A-06 



in the derivation. The independent variables of the dynamic system chosen 
as example satisfy the following conditions: (l) the load is a triangular 
load; (2) the displacement reaches the limiting elastic displacement x 
at time t where ,t < T; and (3) the maximum displacement is reached at 



m 



where 



t < T . 

m 

When t < t , the equation of motion 'is given by equation (5-33a). 

. < t < t . 
e m ' 



When t < t < t , the equation of motion is given by equation (5-33b). 
The expression for 
given by equation (A. 35). 



x derived from these two differential equations is 

m 



m 



m 



+ v 



", t - t 

1 m e 


2 6/T 


- t e)_ 


R 




_a 

2m 


ft - t \ 

\ m e/ 



(A.35) 



where 



v = the velocity of the mass at time t . 

C 6 

This equation shows that x is a function of B, T, 



t and v 

C 6 

pend on B, 



m 
The last three variables: 



m, 



k, .R 



m' 



and 



v , in turn, de- 



m 



k, and R. 



m 



Considering the effect of load duration 



T on x , the general expression of I , 
equation (A. 3), is given by equation (A. 36). 



instead of' being given by 



m 



ox 

m . 

""or 



m 



dx at 

ra m 

. ~5f + 

m 



(A.36) 



After the partial derivatives in equation (6.36) are evaluated, the expres- 



sion of 



X xT iS given by: 




(A. 37) 



By a similar procedure, the expressions for all the percentage -increment 
ratios of an elasto -plastic system subjected to triangular or rectangular 
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loads can be derived. These expressions are summarized in paragraph A-08c. 
The percentage-increment ratios for elasto-plastiu -systems subjected 
to rectangular and triangular loads are given in figures A. 5 through A. 8. 
(The first two f igures are for rectangular loads and the last two, for tri- 
angular loads. In figures A. 5 and A. 7, the percentage-increment ratios 

with respect to T> m and k are plotted against CL with C_ as the 

/ 

running parameter. In figures A. 6 and A. 8, the independent variables are 

B and R - Prom these four figures, the following conclusions can be 



m 



drawn: 



(1) The effect of load duration T . The curves in figures A, 5 and 
A. 7 are divided into two groups by the curve corresponding to CL 1 . 
When C R < 1, the value of I T , with a few exceptions, lies between 1 
and 2. Hence for a given C , the maximum displacement is proportional to 
T at small values of C and is proportional to the square of T at large 
values of C^ . The smaller the values of CL , the smaller will be the 
value of C m when I_ approaches 2. On the other hand when CL > 1.0, 

j; jCU. J\ 

the value of I m lies approximately between zero and unity. Hence, at 
small values of C , the maximum displacement is approximately propor- 
tional to the load duration T , while at large values of C T , the maxi- 
mum displacement is practically independent of the load duration T . 

(2) The effect of the mass m and the spring constant k on the 
maximum displacement can be approximately summarized into the following two 
statements, (l) When C R < 1, the values of 1^ lie approximately be- 
tween -0.5 and -1.0, and the values of 1.^ lie between -0.5 and 0. Hence, 
at small values of CL, , the maximum displacement is inversely proportional 
to the square root of both the mass, m , and the spring constant, k . 

At large values of C T , x ffl is inversely proportional to m but is inde- 
pendent of k . (2) When C R 1.0. At small values of CL, , x ffi is in- 
versely proportional to square root of both m and k . But at large 
values of CL, , x m is inversely proportional to k but is independent 



m 



of m . 



(3) The effect of the peak load B and the plastic resistance R B 
As shown in figures A. 6 and A. 8, the numerical values of 1^ and 1^ 
scatter over a much wider range than the corresponding values of 1^ , 
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X xm r J xk ' F r the convenience of discussion, the curves in these two 



figures can be very roughly divided into three regions, (l) The first 
region covers all cases when C T < 0.6. Practically for all values of CL , 
the values of 1^ lie between 1 and 2, and the values of I He be- 
tween and -1. Thus for low values of C T , x m is directly proportional 
to B but independent of ly while for high values of C , x is pro- 
portional to the square of B and inversely proportional to R . The 
parameters of one story buildings are generally located in this region. 
Since the numerical values of both 1^ and I are not very large, a 
close agreement between the maximum displacement obtained by a design com- 
putation and that of the design specification can generally be achieved. 
(2) The second region corresponds to cases when CL, > 1.0 and C_ > 1.4 or 
C R < 0.4. In this region both 1^ and 1^ are practically constant.. 
The numerical values of neither one is very large. The cases when C < 0.4 
are not of practical importance because the corresponding values of C are 
always greater than 20. On the other hand, the cases when CL > 1.4 gen- 
erally correspond to the region in which the design of structural members 
such as beams, walls or columns are located. (3) The third region cor- 

responds to cases when C_ > 1.0 and 0.4 < CL < 1.4. In this region I 

1 K XB 

is generally greater than 2.5 and I _ is less than -1.5. Therefore, the 

maximum displacement x is very sensitive to any variation in B or R 

m v in 

especially when CL lies between 0.8 and 1.2. In design problems, both B 
n 

and R ffl are subject to uncertainties in their determinations. Therefore 
the value of x obtained by a design computation would not be very re- 
liable. Hence it would be -a good design practice, in general, to select 
the parameter of a design such that the combination of C and C^ lies 
outside of the third region. 

A-0? THE EFFECT OF LOAD SHAPE ON THE RESPONSE OF DYNAMIC SYSTEMS. The 
effect of load shape is partially investigated in EM 1110-345-415 where the 
responses of a dynamic system subjected to different types of load, as shown 
in figure 5.8, are analyzed and presented separately. For each type of load 
with a given geometry, i.e., a triangle or a rectangle, its "shape" or "pro- 
portion" is completely specified by two variables. These are the peak load 
B and the load duration T . When either B or T is varying, the 
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response is also varying continuously. Therefore, in paragraph A-02, pax- 

tial derivatives of x with respect to B or T can be evaluated and 

m 

the effect of B and T can be expressed quantitatively in terms of 
percentage -increment ratios. Here again the percentage-increment ratios for 
rectangular and triangular loads are considered and presented separately. 

If the concept of percentage-increment ratio is to be extended to the 
analysis of the effect of load shape on the response, it would not be suffi- 
cient to consider Just discrete types of load, for example rectangular and 
triangular loads. Instead, another variable has to be added to the load 
pulse so as to form a family of load curves, and the geometrical shape -of 
the load curve will then be varying in a continuous fashion. For example, 
in figure A. 9, when T, varies from to T , the geometrical load shape 

varies continuously from a triangular 
load to a rectangular load. A 
percentage-increment ratio I _ 
could then be defined which expresses 
the relative magnitude of the per- 
centage changes in x and T, 
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Figure A.9. A set of load curves with T} 
as a parameter 

metrical shape of the load on the 
response could then be analyzed quantitatively. 

It is clear that a rigorous investigation of the effect of the load 
shape cannot be carried out without introducing another variable. For an 
elasto-plastic system, there would be 4 non-dimensional variables involved 
and the result could not be conveniently presented in graphical form. Al- 
though for linearly elastic or completely plastic systems there would be no 
difficulty in presenting I _ for a given family of load shapes. Its 
practical application would be limited to the case where the actual load 
shape variation follows the same pattern as the variation in the given 
family. Because of these difficulties, as far as the effect of load shape 
is concerned, a qualitative discussion, supplemented by non-dimensional 
graphs whenever convenient, is presented in this paragraph. Specifically, 
two cases are considered. In the first case, the response of an elasto- 
plastic system to rectangular and triangular loads is compared with each 
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other. This will indicate in what region, if any, the load shape is not 
important. In the second case, the percentage-increment ratio I , given 
by equation (A. 28), for rectangular and triangular loads is presented in 
graphical form. These curves will show how the response is affected by 
changing the "shape" or "proportion" of the load curve without changing its 
geometry. 

a ' Corcpar 3 - 80 " of the Responses of an Elasto-Plastic System Subjected 
to Rectangular and Triangular Loads. In figure A. 10, the ratios of the 
maximum displacement for an elasto-plastic system subjected to a rectangular 
load to the corresponding value for a triangular load are plotted. There 
are three regions in this figure where the shape of the load pulse has no 
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Figure A.10. Ratio of maximum displacements for rectangular to the corresponding values 
in triangular load, elasto-plastic system 
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important effect on the response, (l) In the first region, vhere C T > 0.5 
and C < OA, the ratio of x^s is approximately equal to U. The cor- 
responding values of C and % are always greater than unity. In this 

Jt *- 

region, the maximum strain energy is approximately proportional to x ffi but 
the maximum work done is approximately proportional to the square of the 
area under the load -time curve. Therefore, the maximum displacement is ap. 
proximately proportional to the square of the area under the load curve 
while the detail geometry of the load curve is not important. (2) In the 
second region where C T < 0.5 and C R > 1.2, the ratio of x m 's is approxi- 
mately equal to 2. In this region, x ffl < 2x g , and t ffl > T . Although W B 
is still approximately proportional to" the square of the area of the load 
pulse, the maximum strain energy is also proportional to the square of x . 
There the maximum displacement is proportional to the area under the load- 
time curve and is not much affected by the detail shape of the load. (3) 

In the third region where C m > 2 and CL. > 1.2, the ratio of x 's is ap- 

J. t\ m 

proximately equal to unity. This corresponds to the cases where t is 
small in comparison with T . The maximum displacement is proportional to 
the peak load B while the load shape is not important. 

Outside of the above three regions, the ratios of x 's for rectangu- 
lar and triangular loads vary considerably with both C and C R . No 
simple conclusion can be drawn. Generally the response is expected to be 
affected by the detail geometry or the load curve. 

b. The Percentage-Increment Ratio I T . The percentage-increment 

xo 

ratio 1.^ given by equation (A. 28) is plotted in figures A. 11 and A. 12 
for rectangular and triangular loads respectively. As discussed in para- 
graph A-03c, the value of 1^ indicates the effect of varying the "shape" 
or "proportion" of a load curve without changing its basic geometry. In 

figures A. 11 and A. 12, the numerical values of' I _ are small when 

xj 

C T < 0.3. In this region the response will not be much affected when one 
triangular load is replaced by another triangular load, or one rectangle 
by another rectangle, or any load curve which resembles & triangle by the 
corresponding triangular load. 

A-08 EXPRESSIONS FOR MAXIMUM DISPLACEMENTS AND PERCENTAGE-INCREMENT RATIOS.. 
The expressions of t ffi , x ffl , 1^ and I XT for linearly elastic, 
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Figure A. 11. The percentage~increment ratio l x for an elasto-plaitic 
system subjected to rectangular loads 
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completely plastic and elasto -plastic systems subjected to rectangular and 
triangular loads are given in this paragraph. From the expressions of 
these four 'quantities, the expressions of many other quantities can be de- 
rived. For example, in terms of x ffl , the expression for W ffi is: 

when x < x 
m e 



and when x > x 

m e 



In terms of I , and I , the expressions of many other percentage- 

X J. X\ 

increment ratios can be obtained. These relationships are given by equa- 
tions (A.UO) and (A. 1*1). 



1 - 2C 2C - 1 
x ... _ x 

Y 



1 - 2C 2C - 1 , 

I = I . 1 = -I * T - * f-r ^ 

2 xR xB 2C \m ~ 2C 

X X 

a. Expressions for Linearly Elastic Systems. 
Rectangular Loads 

(i) For t ffl < T 



k x 
D.L.F. = -g-S = 2 

where 

D.L.F. = dynamic load factor 
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(ii) For t > T 



Triangular Loads 

(i) For t < T 



m 



5= - i tan' 1 (a, c 
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(ii) For t m > T 
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b. Expressions for Completely Plastic Systems. 
Rectangular Load 

t , 
ml , , 

m I A k7 N 

rp ft \"^|<l/ 

x 1 - CL 

SL_ .. n R (j 

r^7 C R ^ 

X xT = 2 (' 



2 m 



Triangular Load 



(i) For C_ > 0.5 and t < T 

JK In 

t 

^ = 2 



2m 



(ii) For C B < 0.5 and t > T 

K ffl 

t 

m 
T- 






*\ " 5) 



I xT - 2 (A.49c) 

^TT^ 



I - 2 (A.lf9c) 
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c. Exrnressions for Elasto-Plastic Systems. 
Rectangular Load 

(i) For x ffl < x g and t ffl < T . .The expressions for this case are 
the same as those given by equation (A.lj-2). 

(ii) For x ffl < x g and t ffl > T . The expressions for this case are 
the same as those given by equation (A.^3). 

(iii) For x > x, t > T and t > T . 
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- cos 
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(A.50a) 
(A.50b) 
(A.50c) 
(A. 501) 
(A.50e) 



(iv) For x ffl > x e , t e < T and t ffl < T . 
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(A.51a) 
(A.5Ib) 
(A.51c) 
(A. 514) 
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(v) For x m > x e , t e < T and t ffl > T . 

t 



T = Sr 
n 



(A.52a) 
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Triangular Load 

(i) ^or x < x and t < T . The expressions are the same as 
x ' m e ro 

those given by equation ( 

(ii.) For x < x , and t > T . Tlie expressions are the same as 
^ ' me' m 



those given by eq_uation (A.^5). 

(iii) For x m > X Q , t a > T, a.nd 



T . 



sin 2 ~ m 2 C R 
n 



/ %> 

cos J2n 7=- 

V r>y 



sin 2* 



.) 



C T ! (A.53a) 



** *e 1 

_ 1 __^ J ^^^ j. 

T ~ T - li* 2 r 

:i n ^rt 0. 



cos 2n 7p 4- 2^ T sin 
n 



oos 2rt ~ - C r 



x 
J 
x 



ft 






x \ n n 



xR 



^!L/J2 !e 

c IT " ^ 

x \ n ~n 



(A.53c) 

(A. 53d) 
(A.?3e) 



A-08c 



EM 1110-3^5-415 
15 Mar 57 



(iv) For x m > x e , t < T, and t < T 
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